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^xhe  plane  wave  reflection  coefficient  R  of  a  horizontally  stratified  ocean 
bottom  is  studied  for  typical  deep  sea  sediment  types  in  order  to  determine 
the  importance  of  sediment  rigidity  (shear  wave  propagation)  on  bottom 
reflection  loss,  RL  =  -20  log^Q(|R|)  dB.  The  emphasis  is  on,  but  not 
restricted  to,  low  grazing  angles  (0  to  433’)  and  low  frequencies 
(10  to  200  Hz),  R  is  studied  using  a  computational  model  accurate  for  fre¬ 
quencies  above  10  Hz  which  treats  a  single  inhomogeneous  (solid)  sediment  _ 
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layer  overlying  a  homogeneous  (solid)  substrate ^ Comparing  RI.  due  to  a  solid 
sediment  RLS  with  RL  oi  a  fluid  sediment  RLE  yields  ARL=R1.S-RLF,  the  change  in 
RL  produced  by  sediment  rigidity. 

For  a  hypothetical  turbidite  layer  the  dependence  of  ARL  on  layer  thickness 
suggests  an  empirical  classification  of  sediments  into  "thick"  and  "thin" 
categories.  For  a  thick  (500  m)  layer  ARl.'O  at  low  grazing  angles  and  sediment 
rigidity  can  be  neglected.  For  a  thin  (36  m)  layer  at  a  grazing  angle  of  20°, 

ARL  -  30  dB  at  20  Hz  and  ARL  =  4  dB  at  200  Hz;  sediment  rigidity  in  this  case  is 
not  negligible  but  produces  significant  increases  in  reflection  loss. 

The  frequency  dependence  of  RLS  for  the  36  m  thick  layer  reveals  high  and  low 
frequency  regimes  with  strikingly  different  behavior.  For  a  grazing  angle  of 
20°,  ARL  is  about  2  dB  at  low  frequencies  with  very  large  peaks  (-25  dB) 
occurring  at  2  Hz  intervals.  These  peaks  are  the  result  of  interference  effects 
related  to  sediment  shear  wave  propagation.  The  interference  peaks  decrease  in 
magnitude  at  higher  frequencies  as  shear  wave  attenuation  increases  and  reduces 
the  magnitude  of  interference  effects.  At  high  frequencies  the  peaks  have 
entirely  disappeared.  A  value  of  ARL  ~  4  dB  persists  and  there  is  also  a  gentle 
oscillation  occurring  in  both  RLS  and  RLF  due  to  sediment  compressional  wave 
interference  effects. 

An  examination  of  the  reflection  and  transmission  coefficients  at  the  water- 
sediment  and  sediment-substrate  interfaces,  by  means  of  an  expansion  in  the  ratio 
of  sediment  shear  speed  to  the  sound  speed  in  water,  reveals  the  mechanism  for 
exciting  sediment  shear  waves  and  provides  an  understanding  of  the  empirical 
separation  into  ''thick''  and  ''thin*'  sediments.  At  the  water-sediment  interface 
sediment  shear  waves  are  not  efficiently  excited  and  the  sediment  can  be  treated, 
in  most  cases,  as  a  fluid.  The  dominant  mechanism  for  exciting  sediment  shear 
waves  is  compressional  wave  conversion  at  the  sediment-substrate  interface.  At 
a  given  grazing  angle  in  a  thick  sediment  the  compressional  wave  will  have  a 
turning  point  well  above  the  sediment-substrate.  Compressional  wave  amplitude 
at  the  Interface  will  be  small.  As  a  result  very  little  energy  is  transferred 
[to  sediment  shear  waves,  and  th^\sediment  can  be  accurately  described  as  a  fluid. 
|For  a  thin  sediment  the  compressional  wave  strikes  the  interface  or  has  its 
turning  point  near  the  interface;  cotnpressional  wave  amplitude  at  the  interface 
is  significant,  sediment  shear  waves  are  excited,  and  the  sediment  must  be 
treated  as  a  solid.  The  difference  between  a  "thick"  and  a  "thin"  layer  is  then 
translated  into  a  question  of  whether  the  compressional  wave  significantly 
interacts  (thin)  with  the  substrate  or  not  (thick). 


Parameter  studies  are  reported  which  investigate  the  dependence  of  RLS  on  bottom 
parameters  and  gradients.  Homogeneous  sediment  layers  and  layers  with  gradients 
are  studied  to  determine  which  parameters  are  important  in  describing  the  acoustic 
properties  of  solid  sediments.  Variations  about  a  parameter  set  indicate  the 
accuracy  needed  in  specifying  bottom  parameters  to  accurately  predict  reflection 
loss . 
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I.  INTRODUCTION 


Acoustic  interaction  with  the  ocean  bottom  is  an  important  factor  in 
determining  the  propagation  of  low  frequency  sound  waves  in  the  ocean. 
Because  of  the  long  wavelengths  involved,  the  acoustic  energy  penetrates 
into  the  ocean  bottom  and  interacts  with  the  subbottom  structure.  The 
excitation  of  compressional  (P),  shear  (S),  and  interface  waves  in  the 
sediment-substrate  structure  are  important  energy  loss  processes 
affecting  the  propagation  of  low  frequency  sound,  particularly  when 
compared  to  the  negligible  attenuation  occurring  within  the  water  column 
itself. 

Determining  the  dependence  of  the  plane  wave  reflection  coefficient 

R  (and  its  analog,  the  reflection  loss  RL  =  -20  log^^ ( | R | ) )  on  subbottom 

parameters  is  a  useful  and  traditional  method  of  studying  the  bottom 

interaction.  Parameter  studies  of  R,  using  computational  models  of  the 

bottom  interaction  for  an  assumed  fluid  sediment,  have  determined  the 

importance  of  several  subbottom  features.  Density  gradients  have  been 

found  to  be  negligible  in  typical  deep  sea  sediments.'*'  The  effect  of 
2  3  4 

P  wave  speed  ’  and  absorption  profiles  have  been  treated.  Interface 
wave  excitation  at  the  sediment-substrate  interface  has  been  identified 
as  a  loss  mechanism. The  effect  of  surface  roughness  at  the  water- 
sediment  interface  has  also  been  studied. ^ 

The  effect  of  sediment  rigidity  (shear  wave  propagation)  on  R  has 

received  little  attention  and  is  not  well  understood.  Although  included 

8  9 

in  some  computational  models,  ’  no  systematic  studies  have  been  made  to 

determine  the  importance  of  sediment  S  wave  propagation.  A  detailed  study 

at  this  time  would  provide  realistic  and  useful  results  since  parameter 

ranges  can  now  be  constrained  by  recently  available  estimates  of  sediment 

S  wave  properties. Such  a  study  is  of  importance  since  the  relatively 

high  attenuation  of  sediment  S  waves  makes  their  excitation  a  potentially 

13  14 

important  energy  loss  mechanism.  Recent  work  *  indicates  that  sediment 
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S  wave  excitation,  while  negligible  in  thick  sediment  layers,  can  be 
the  dominant  energy  loss  mechanism  in  thin  sediment  layers. 

In  this  paper  we  present  the  results  of  our  investigation  of  the 
effect  of  sediment  rigidity  on  R.  Our  emphasis  is  on,  but  not  restricted 
to,  low  grazing  angles  and  typical  deep  sea  sediment  types.  Our  study 
shows  that  the  mechanism  for  sediment  S  wave  excitation  is  primarily  P 
wave  conversion  at  the  sediment-substrate  interface.  (The  gradient- 
driven  P-S  conversion  is  negligible  for  these  sediment  types  for  fre¬ 
quencies  above  about  3  Hz.^)  An  examination  of  the  reflection  and 
transmission  coefficients  shows  that  this  is  due  to  the  small  surficial 
shear  wave  speed  in  the  sediment.  Typical  deep  sea  sediments  can  be 
accurately  treated  as  fluids  at  the  water-sediment  interface  where 
sediment  S  wave  excitation  is  negligible.  This  excitation  mechanism  also 
allows  a  classification  of  sediment  layers  into  "thick"  and  "thin" 
categories.  For  a  thick  sediment  the  P  wave  has  its  turning  point  well 
above  the  sediment-substrate  interface.  Consequently,  sediment  S  wave 
excitation  is  negligible  and  the  sediment  can  be  accurately  treated  as  a 
fluid.  For  a  thin  sediment  the  P  wave  significantly  interacts  with  the 
substrate.  Shear  waves  are  excited  in  the  sediment,  and  the  sediment 
must  be  treated  as  a  solid.  For  thin  sediments  our  study  shows  low  and 
high  frequency  regimes  with  strikingly  different  behavior.  At  low 
frequencies  sediment  shear  wave  attenuation  is  small,  and  interference 
effects  due  to  propagation  of  S  waves  within  the  sediment  column  dominate 
the  dependence  of  R  on  frequency.  Large  (>20  dB)  peaks  in  RL  occur.  The 
frequency  separation  between  these  peaks  is  related  to  changes  in  S  wave 
phase  of  2n.  At  higher  frequencies  the  increased  shear  wave  attenuation 
destroys  this  peak  structure.  The  shear  waves  excited  at  the  sediment- 
substrate  interface  are  then  totally  absorbed  within  the  sediment,  result¬ 
ing  in  an  almost  constant  additional  loss  (-4  dB)  compared  to  the  RL  of  a 
fluid  sediment.  We  also  present  parameter  studies  aimed  at  identifying 
important  parameters  in  solid  sediments  and  the  accuracy  to  which  these 
parameters  must  be  known  to  accurately  compute  RL. 
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The  cool  used  in  this  study  is  a  recently  developed  computational 
model1^  of  R  which  allows  sediment  shear  wave  propagation.  The  model 
assumes  horizontal  stratification  and  treats  a  single  sediment  layer 
overlying  a  semiinfinite  homogeneous  substrate.  The  model  is  based  on 
numerical  integration  of  the  depth  separated  wave  equations  for  the 
potentials  giving  rise  to  the  compressional  and  shear  waves  in  the 
sediment,  and  thus  includes  all  wave  properties,  sue)  is  penetration 
beyond  turning  points,  boundary  wave  excitation,  et  The  use  of 
numerical  integration  permits  sediment  parameters  to  have  arbitrary 
depth  dependencies  and  to  be  individually  varied  to  determine  their 
influence  on  R.  The  basic  approximation  of  the  model  is  the  use  of 
the  Helmholtz  equations  with  depth  dependent  wave  number  to  describe 
the  potentials.  This  is  essentially  a  high  frequency  approximation 
in  which  effects  directly  dependent  on  gradients  and  the  continuous 
coupling  between  shear  and  compressional  waves  are  neglected.  For 
parameters  typical  of  deep  sea  sediments,  this  model  is  accurate  for 
frequencies  above  about  10  Hz. 

The  first  section  of  this  paper  discusses  a  hypothetical  turbidite 
layer  with  realistic  parameters  and  depth  variations.  By  varying  the 
sediment  layer  thickness  and  the  wave  frequency,  our  basic  results  are 
observed.  The  second  section  develops  a  theory  of  the  reflection  and 
transmission  coefficients  based  on  an  expansion  in  the  sediment  S  wave 
velocity.  These  are  used  to  construct  a  ray  picture  of  the  acoustic 
interaction,  which  is  then  used  to  interpret  the  results  of  the  hypo¬ 
thetical  turbidite  layer  study.  The  wave  theory  model  verifies  some 
parameter  dependencies  predicted  by  the  ray  model.  The  next  section 
presents  the  results  of  parameter  studies  for  homogeneous  layers  and 
layers  with  gradients.  The  parameters  of  a  realistic  layer  are 
perturbed  to  obtain  an  estimate  of  the  accuracy  required  in  the  values 
of  these  parameters  to  compute  R  accurately.  The  effect  of  sediment 
rigidity  on  interface  wave  excitation  is  also  treated.  The  final 
section  summarizes  the  results  of  our  study. 
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II.  HYPOTHETICAL  TURBIDITE  LAYER 


In  this  section  we  consider  the  effect  of  sediment  rigidity  on  the 

reflection  less  from  a  hypothetical  turbidite  layer.  RL  is  studied  as 

a  function  of  grazing  angle  6,  frequency  f,  and  sediment  layer  thickness 

H  (the  grazing  angle  is  the  complement  of  the  angle  of  incidence  , ,  i.e., 

9 

-’  =  90-r).  Our  results  agree  with  the  recent  work  of  Fryer  for  a  thick 
sediment  layer  and  illus' rate  our  most  important  results  for  thin 
sediment  layers. 

The  depth  structure  of  the  turbidite  layer  is  given  in  Table  I. 

9 

These  parameters  were  obtained  from  the  recent  work  of  Fryer.  Constant 
gradients  are  assumed  between  the  depths  given  in  Table  I,  the  parameter 
values  are  representative  of  this  type  of  sediment,  and  the  attenuation 
is  assumed  to  depend  linearly  on  frequency.  Fryer's  650  m  sediment 
layer  was  truncated  at  518  m  so  that  our  numerical  model  would  produce 
accurate  results  at  low  grazing  angles  for  a  frequency  of  20  Hz.  For  a 
significantly  thicker  layer  (650  m)  the  S  wave  attenuation  and  the 
exponential  behavior  of  the  P  wave  potential  combine  with  the  finite 
computer  word  length  to  produce  a  loss  of  precision  for  grazing  angles 
below  about  20°.  The  accuracy  of  the  calculated  RL  at  518  m  was 
established  in  two  ways.  First,  for  a  solid  sediment,  RL  below  9=25° 
agrees  to  within  0.1%  for  H  =  518  m,  385  m,  and  252  m.  Second,  for  the 
same  thicknesses  and  grazing  angles,  RL  for  both  fluid  and  solid 
sediment  structure  agreed  to  within  1%.  The  theoretical  analysis  of 
the  next  section  of  this  report  shows  that  this  constancy  of  RL  with 
depth  at  low  grazing  angles  and  the  lack  of  dependence  on  sediment 
rigidity  are  expected  for  these  relatively  thick  sediments.  Since  it 
is  less  and  less  difficult  to  maintain  precision  as  H  decreases,  the 
accuracy  of  the  computation  at  518  m  is  verified. 


TABLE  I 


ACOUSTIC  PARAMETERS  OF  THE  HYPOTHETICAL  TURBIDITE  LAYER 

c  is  the  compressional  wave  speed; 

v  is  the  shear  wave  speed; 

p  is  the  density; 

kp  is  the  compressional  wave  attenuation; 

kg  is  the  shear  wave  attenuation. 


Depth 

(m) 

c 

(m/sec) 

k 

P 

(dB/m/kHz) 

P 

(g/cm3) 

V 

(m/sec) 

k 

s 

(dB/m/kHz) 

Water 

1530 

— 

1.030 

i 

1 

1 

0 

0.065 

1.530 

116 

8.46 

36 

0.100 

1.579 

283 

5.60 

120 

1674 

0.200 

1.689 

391 

8.60 

518 

1992 

0.135 

2.010 

621 

4.35 

Substrate 

4460 

0.008 

2.460 

2400 

0.04 
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Figure  1  shows  RL  as  a  function  of  9  for  f  =  20  Hz  and  H  =  518  m. 

The  three  curves  were  obtained  for  different  sediment-substrate  config¬ 
urations.  The  dotted  line  was  obtained  for  a  fluid  sediment  and  fluid 
substrate  (the  FF  case),  the  dashed  line  for  a  fluid  sediment  and  solid 
substrate  (the  FS  case),  and  the  solid  line  for  a  solid  sediment  and 
solid  substrate  (the  SS  case).  The  remaining  case  of  a  solid  sediment 
and  a  fluid  substrate  (the  SF  case)  is  not  shown  in  Fig.  1.  It  is 
essentially  identical  to  the  FF  case  over  the  entire  range  of  grazing 
angles,  frequencies,  and  sediment  thicknesses  considered.  The  parameters 
for  the  fluid  sediment  and  fluid  substrate  were  obtained  from  Table  I 
by  ignoring  the  S  wave  parameters. 

Comparing  the  curves  of  Fig.  1  shows  the  effect  of  sediment  rigidity 
in  this  thick  layer.  Case  FF  has  the  expected  small  RL  for  small  9 
followed  by  an  oscillatory  structure  due  to  interference  between  the  wave 
directly  reflected  from  the  water-sediment  interface  and  the  P  wave  pene¬ 
trating  the  sediment  and  returning  to  the  water.  Beyond  the  critical 

angle  for  substrate  P  wave  propagation  (0  =70°),  RL  increases  as  energy 

P 

is  carried  away  by  substrate  P  waves.  Comparing  the  FS  and  FF  cases 

shows  the  effect  of  substrate  rigidity.  At  both  small  and  large  grazing 

angles  the  FS  and  FF  cases  are  nearly  identical,  an  indication  that 

substrate  rigidity  has  a  negligible  effect  at  these  angles.  RL,  however, 

increases  between  the  critical  angle  for  substrate  S  wave  propagation 

0g  and  0p.  Substrate  S  waves  carry  away  additional  energy  in  this  angular 

range.  Comparing  the  SS  and  FS  cases  shows  the  effect  of  sediment 

rigidity  alone.  Except  for  grazing  angles  between  0  and  0  ,  sediment  S 

s  p 

wave  propagation  has  a  negligible  effect  on  RL.  The  major  effect  between 
the  critical  angles  is  a  shift  in  the  peak  structure. 

Comparison  of  Fig.  1  with  Fryer's  Fig.  8  shows  qualitative  agreement 
between  the  corresponding  FF  and  SS  cases.  Fryer  obtained  his  result 
using  a  different  computational  technique  and  for  a  thicker  (650  m) 
sediment  layer.  The  conclusion  can  be  drawn  from  our  intermediate  FS 
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that  sediment  rigidity  is  negligible  for  this  particular  sediment 
structure  at  low  grazing  angles  and  at  a  frequency  of  20  Hz. 

Figure  2  shows  RL  for  a  thinner  sediment  layer,  H  =  252  m,  at 

f  =  20  Hz.  The  FS  and  SS  cases  are  shown.  Comparison  of  the  solid  (SS) 

and  dashed  (FS)  curves  shows  that  the  presence  of  sediment  rigidity 

produces  the  expected  shift  in  the  oscillatory  structure  between  the 

critical  angles,  0  and  0  .  In  addition,  the  SS  case  has  an  increased 
s  p 

RL  beginning  at  a  minimum  grazing  angle,  0  =24°,  and  extending  up  to 

o 

about  45°.  This  increase  in  RL  is  entirely  due  to  sediment  rigidity. 

For  grazing  angles  below  0^  the  effect  of  sediment  rigidity  is 
negligible. 

Figure  3  shows  RL  for  the  252  m  layer  as  a  function  of  frequency 
at  0=2O°<0o.  The  frequency  range  is  10  to  100  Hz.  A  comparison  of 
the  SS  and  FS  cases  shows  that  for  frequencies  above  15  Hz  sediment 
rigidity  is  negligible.  Below  15  Hz  the  difference  between  the  curves 
may  not  be  meaningful  since  this  is  near  the  low  frequency  limit  of 
the  validity  of  the  computational  model. 

A  consideration  of  RL  as  a  function  of  0  for  sediment  thicknesses 

less  than  252  m  reveals  a  relationship  between  H  and  0^.  As  H  decreases, 

9  also  decreases.  Thus,  as  sediment  thickness  decreases,  sediment 
o 

rigidity  affects  RL  at  steadily  lower  grazing  angles.  In  fact,  sediment 
S  wave  excitation  can  be  the  dominant  energy  loss  mechanism  in  thin 
sediment  layers.  This  is  illustrated  in  Fig.  4  which  compares  RL  for 
the  SS  and  FS  cases  at  H  =  36  m  and  f  =  20  Hz.  The  angle  0^  is  0°. 
Sediment  S  wave  excitation  is  responsible  for  the  dramatic  20  dB 
increase  in  RL  near  0=12°. 

The  dependence  of  RL  on  frequency  for  the  same  36  ra  layer  is  shown 
in  Fig.  5.  The  SS  and  FS  cases  are  compared  from  10  to  200  Hz  at  0=20°. 
The  SS  case  has  high  and  low  frequency  regimes  with  strikingly  different 
properties.  At  low  frequencies  the  effects  of  sediment  rigidity  increase 
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REFLECTION  LOSS  -  dB 


FIGURE  3 

REFLECTION  LOSS  versus  FREQUENCY  FOR  A  252  m 
THICK  HYPOTHETICAL  TURBIDITE  LAYER  AT  20° 
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the  level  of  RL,  relative  to  the  FS  case,  by  about  3  dB,  with  additional 
very  large  peaks  of  up  to  25  dB.  The  peaks  repeat  with  an  almost  constant 
separation  of  2.65  Hz.  The  magnitude  of  the  peaks  generally  decreases 
with  increasing  frequency;  however,  any  individual  peak  may  be  higher  or 
lower  than  neighboring  peaks.  The  peaks  disappear  by  about  90  Hz  and 
the  high  frequency  regime,  characterized  by  rather  broad  oscillations  in  RL, 
begins.  The  maxima  of  these  broad  oscillations  are  separated  in  frequency 
by  about  65.5  Hz  and  occur  in  both  the  SS  and  FS  cases.  In  this  high 
frequency  regime,  sediment  S  wave  excitation  produces  an  almost  constant 
increase  in  RL  of  about  4  dB. 

The  large  peaks  in  the  low  frequency  regime  of  Fig.  5  are  associated 
with  interference  effects  due  to  sediment  S  wave  propagation.  This 
identification  is  made  by  noting  that  the  frequency  interval  between  peaks 
is  related  to  S  wave  parameters.  One  expects  interference  effects  to 
recur  when  the  relative  phases  of  the  waves  involved  change  by  2n .  The 
phase  of  the  S  wave,  accumulated  in  traveling  through  the  sediment  and 
back,  changes  by  2^  for  a  frequency  change  Af=v2/2H,  where  v2  is  the 
average  sediment  S  wave  speed.  (The  sediment  shear  speed  is  so  small 
that  the  S  waves  propagate  nearly  normal  to  the  plane  of  stratification.) 

The  measured  separation  between  peaks  is  within  5%  of  the  estimated 
value  of  Af.  The  corresponding  quantity  for  P  waves  (including  the 
grazing  angle  correction)  estimates  a  frequency  change  of  65.4  Hz,  which 
is  p  factor  of  30  too  large  to  explain  the  peak  separation.  In  addition, 
the  decrease  in  peak  height  with  frequency  appears  to  be  related  to  the 
increased  attenuation  of  S  waves  with  frequency.  At  90  Hz  the  attenua¬ 
tion  is  estimated  to  be  23  dB  for  an  S  wave  traveling  one  way  through 
the  sediment.  This  calculation  predicts  that  S  wave  interference  would 
be  negligible  above  90  Hz  (maximum  effect  of  1%). 

The  gentle  oscillations  above  90  Hz  in  Fig.  5  are  due  to  P  wave 


interference  effects.  This  conclusion  is  reached  because  the  same 
oscillations  are  observed  in  both  the  SS  and  FS  cases.  The  measured 


separation  between  maxima  is  65.5  Hz,  which  agrees  very  well  with  the 
estimated  Af  =  65.4  Hz. 

Figure  6  shows  that  sediment  rigidity  is  still  significant  in  the 

high  frequency  regime.  The  SS  and  FS  cases  are  compared  at  H  =  36  m 

and  f  =  200  Hz.  The  increase  in  RL  due  to  sediment  S  wave  excitation  is 

significant  (-4  dB)  between  15°  and  30°.  The  minimum  grazing  angle, 

9  ,  is  about  5°. 
o’ 

Inspection  of  Fig.  5  near  20  Hz  shows  that  the  very  large  increase 
in  RL  seen  in  the  SS  case  in  Fig.  4  is  due  to  the  fortuitous  location 
of  one  of  the  interference  peaks  at  20  Hz.  For  comparison  with  Fig.  4, 

RL  is  given  at  21  Hz  (minimum  in  RL)  in  Fig.  7.  The  increase  in  RL 
due  to  sediment  S  wave  excitation  is  still  present  at  low  angles  but  its 
magnitude  is  2  dB  rather  than  20  dB.  This  is  still  a  significant  increase 
in  RL  when  compared  to  the  RL=0  for  the  fluid  sedimei.t  case. 

The  dependence  of  6^  on  H  gives  an  indication  of  the  physical 

mechanism  responsible  for  the  increase  in  RL  due  to  sediment  rigidity. 

For  a  given  grazing  angle  the  compressional  wave  in  the  sediment  has  a 

turning  point  at  a  depth  H^.  For  a  thick  sediment  layer  H^  will  be  well 

above  the  sediment-substrate  interface.  As  the  grazing  angle  increases, 

H  increases.  At  some  critical  angle  0  ,  H  =H  and  the  compressional  wave 
t  c  t 

interacts  with  the  sediment-substrate  interface.  As  the  thickness  of  the 

sediment  layer  decreases,  the  angle  9  also  decreases.  This  is  the  same 

c 

qualitative  behavior  exhibited  by  the  angle  9q,  above  which  sediment 
rigidity  causes  an  increase  in  RL.  This  similarity  suggests  that  sediment 
shear  wave  excitation  by  compressional  wave  conversion  at  the  sediment- 
substrate  interface  produces  the  increase  in  RL. 

Two  other  observations  support  this  identification  of  P  wave 
conversion  at  the  sediment-substrate  interface  as  the  dominant  mechanism 
for  sediment  S  wave  excitation.  One  concerns  the  negligible  influence  of 
the  water-sediment  interface  while  the  other  shows  the  importance  of  the 
substrate  interface.  S  wave  excitation  at  the  water-sediment  interface 
appears  to  be  negligible.  The  properties  of  the  water  interface  remained 
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constant  while  the  sediment  thickness  varied  and  large  changes  in  RL 
due  to  sediment  rigidity  occurred.  It  is  difficult  to  establish  a  causal 
relationship  between  something  which  remains  constant  and  another  quantit 
which  changes  significantly.  The  importance  of  the  substrate  interface  i 
further  indicated  by  the  negligible  difference  between  RL  for  the  FF  and 
SF  cases  and  the  large  changes  seen  between  the  SS  and  FS  cases.  This 
indicates  that  the  nature  of  the  substrate  is  important.  Since  the 
substrate  enters  only  through  the  sediment-substrate  interface  conditions 
this  indicates  that  sediment  S  waves  originate  at  the  substrate  interface 

The  results  of  this  section  lead  to  an  empirical  classification  of 
sediments  into  thick  and  thin  lavers  based  on  the  effect  of  sediment 
rigidity.  The  behavior  of  thin  layers  can  be  further  split  into  high 
frequency  and  low  frequency  regimes,  based  on  S  wave  attenuation.  For 
thick  sediments  little  energy  is  transferred  to  sediment  S  waves  and  the 
sediment  can  accurately  be  treated  as  a  fluid  at  low  grazing  angles.  For 
a  thin  sediment  significant  energy  can  appear  in  sediment  b  waves.  Sedi¬ 
ment  rigidity  is  not  negligible  and  the  sediment  must  be  treated  as  a 
solid.  For  thin  sediments  at  low  frequencies  (S  wave  attenuation  less 
than  20  dB  across  the  layer)  interference  effects  due  to  S  wave  propa¬ 
gation  are  important.  Large  (>20  dB)  increases  in  RL  are  possible  at 
particular  frequencies.  At  high  frequencies  (S  wave  attenuation  larger 
than  20  dB  across  the  layer)  sediment  S  wave  propagation  can  be  neglected 
Significant  energy  can  still  be  lost  to  S  waves  that  are  totally 
absorbed  within  the  sediment. 
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III.  MECHANISM  OF  SEDIMENT  SHEAR  WAVE  EXCITATION 


In  this  section  a  qualitative  ray  model  of  P  and  S  wave  propagation 
in  a  solid  sediment  layer  is  developed.  The  qualitative  predictions  of 
this  model  are  tested  using  the  wave  theory  computational  model.  The 
classification  scheme  of  the  previous  section  is  predicted  by  the  ray 
model  and  is  made  more  understandable.  Basic  to  the  construction  of 
the  ray  model  is  an  understanding  of  the  effect  of  sediment  rigidity  on 
the  reflection  and  transmission  coefficients  at  the  water-sediment  and 
sediment-substrate  interfaces.  These  coefficients  are  examined  using 
a  small  parameter  expansion  based  on  the  relatively  small  S  wave  speeds 
in  marine  sediments.  This  expansion  theory  leads  to  a  detailed  under¬ 
standing  of  the  excitation  mechanism  of  sediment  S  waves. 

A.  Expansion  Theory  of  Reflection  and  Transmission  Coefficients 


The  following  conventions  and  definitions  are  used.  Subscripts  1, 

2,  and  3  identify  quantities  in  the  water,  sediment,  and  substrate, 

respectively.  For  simplicity  we  assume  a  horizontally  stratified 

ocean  bottom  with  homogeneous  media.  A  harmonic  time  dependence  with 

angular  frequency  u  is  assumed.  Compressional  wave  speeds  are  c.,  S 

wave  speeds  are  v.,  and  densities  are  p..  Wave  numbers  are 

J  j  2  21/2  9  9  1 / 9 

defined  as:  K.=w/c.,  B.=w/v.,  k=Kcos8,  <  =[K-k  ]  ,  and  8  . = [ B . -k  ] 

J  J  J  J  JJ  JJ 

where  0  is  the  grazing  angle  of  the  incident  ray  in  the  water.  Subscripts 

p  and  s  will  identify  reflection  and  transmission  coefficients  for  P  and 

S  waves,  respectively.  The  water-sediment  interface  will  be  indicated 

by  W/S  and  the  sediment-substrate  interface  by  S/S. 


The  z  axis  increases  downward  from  the  water  into  the  sediment.  The 
displacement  u.  is  defined  in  terms  of  a  scalar  potential,  ■?  .=<t> .  (z)e^ ^  X 


and  a  vector  potential,  ^=yXj(z)e 


i(kx-wt) 


,  by  u  ,=V4> .+  VxY  . ,  where 
-J  -  J - J 
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F 


y  is  a  unit  vector  along  the  y  axis.  'Hu*  depth  dependencies  are  given 
.  -Ik.z  iliz 


+  _lk  tz 

by  ^-w'e  1  +wje  J  ,  and  X  -V^e 


^  | Z ,  where  W,  and  V 
+V.e  |  J 

J  +  f 

are  the  amplitudes  of  the  downward  propagating  waves  and  W^  and  are 

the  amplitudes  of  the  upgoing  plane  waves. 


The  basis  of  our  study  ot  the  reflection  and  transmission 
coefficients  is  an  expansion  in  the  small  parameter  t  ”v.,/c  .  Only  the 
lowest  order  correction  terms  in  i  will  be  kept.  This  will  be  an 

excellent  approximation  at  the  W/S  Interface  where  the  surflclal  shear 

wave  velocities  are  between  50  m/ sec  and  150  tu/sec,  giving  O.ObtMl.lO. 
The  expansion  parameter  is  not  as  good  at  the  S/S  Interface  because  ot 
the  larger  S  wave  velocity  due  to  positive  shear  speed  gradients.  ^ 
Higher  order  corrections  would  then  he  needed  for  accurate  values  ot 
the  reflection  and  transmission  coefficients  at  the  S/S  interface. 
However,  the  lowest  order  corrections  contain  the  qualitative  results  we 
are  seeking  and  are  adequate  for  thin  layers.  For  the  turbidite  layer 

of  the  previous  section  i  varies  from  0.076  at  the  W/S  Interface  to 

0.185  at  16  m,  0.256  at  120  m,  and  0.406  at  518  m.  liven  at  518  m  the 
second  order  terms  contribute  only  15*  corrections  compared  to  the  40* 
corrections  due  to  first  order  terms. 


Beginning  at  the  W/S  interface  with  a  unity  amplitude  sound  wave 

ik  j  7  -  ik  |  7.  kx  v 

incident  from  the  water,  ^-(o  +R  e  )e  '  *  we  obtain 


K  -  ■)  V°M°>  2 

UP 1 2  R12  ~  2'  1 12  *21  008  ' 


(1) 


T  .  4-  )  Ki0)t0)  2 

%12  T12  +  2‘  12  R21  108  0 


(2) 


Ts12  “  2,2'rj2)(k2/k)  ‘'0s2° 


O) 


where  r|^-(k  -  xyi  ^)/ (KjP  ^+KjP  t )  and  '1^- 2x  ^  J  (k  j+ky>  \  >  aro  lho 

zero  order  reflection  and  transmission  coefficients  obtained  tor  a  iluld 
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sediment  with  the  P  wave  incident  from  medium  i. 


Equations  (1)  through  (3)  show  that,  as  far  as  the  incident  sound 

wave  is  concerned,  the  sediment  can  be  accurately  treated  as  a  fluid. 

Equations  (1)  and  (2)  show  that  sediment  rigidity  produces  only  second 
2 

order  (e  )  corrections  to  the  fluid-fluid  interface  reflection  and 

transmission  coefficients.  Equation  (3)  shows  that  the  excited  sediment 

2 

shear  wave  amplitude  is  of  order  e  and  thus  contains  little  incident 
energy. 

The  downward  propagating  P  and  S  waves  next  strike  the  S/S  interface. 
For  the  incident  P  wave  we  find 


R  ..  *  -  eT  n  cos0 

p23  p23  p 


-  eT  t  cose 
p23  p23  p  p 


R  ®  -2eF  cos0  , 
s23  p 


-  err  cose 
s23  s23  p  s 


For  an  incident  S  wave  we  have 


rP23  =  nd-erscose) 


tp23  =  TP(1~ersCOS0) 


r  „.  =  1  -  2eT  cose  , 
s23  s 


Cs23  =  Ts(1~ersCOse) 


In  Eqs.  4  through  (11)  we  have  used  the  quantities 


(Rp23-Tp23+63T»23/k)/2 


(4) 

(5) 

(6) 

(7) 


(8) 

(9) 

(10) 

(11) 


(12) 
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(13) 


r  =  r  +1/2 

P  s 

n  =  2kJ  /k 
P  2 


[(l-R<°>)p  ■ 

/o  K  T(°>  +  T*-0^  6  /  k 

|\  p23A2 

2/3  3  p23  s23  2  3/3 

*  (1-Rp23) 


d2^P3  +  Ts23k/P2  +  Tp23 


(U) 

(13) 

(16) 


and  the  zero  order  reflection  and  transmission  coefficients  for  a  P  wave, 
with  horizontal  wave  number  k,  incident  from  a  fluid  sediment^ * 


T(o) 

P2  3 


2P2K2^1-2a  ^ 


A 


(17) 


(18) 


,  2 
=  4p2<2a  k 


(19) 


where  a=v.  cosO/c  ,  and 
3  o 


D  =  P3<2 


j^l-2a2j  +  4a4<203y/k2  ]  +  pn< 
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(20) 


Equations  (4)  through  (11)  show  that  sediment  rigidity  can  be 

important  at  the  S/S  interface.  According  to  Eq.  (6)  the  P  wave 

incident  upon  the  S/S  interface  excites  a  first  order  (e^)  S  wave  in 

the  sediment.  Equations  (4),  (5),  and  (7)  show  that  the  reflected  P 

wave  and  P  and  S  waves  transmitted  into  the  substrate  also  have  first 

order  corrections  to  their  amplitudes.  Recalling  that  the  S  wave 

2 

excited  at  the  W/S  interface  is  of  order  e  [Eq.  (3)],  Eqs.  (8)  through 

(10)  show  that  the  reflected  and  transmitted  waves  due  to  it  are  at 

2  1 
most  of  order  e  .  These  are  then  a  correction  to  the  larger,  order  e  , 

corrections  due  to  P  wave  conversion  at  the  S/S  interface.  The  presence 

of  gradients  will  increase  the  relative  effectiveness  of  the  S/S  interface. 
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Equations  (8)  through  (11)  show  an  apparent  inconsistency.  The 
incident  S  wave  is  perfectly  reflected  to  order  e°  and  yet  sets  up 
other  reflected  and  transmitted  waves  whose  amplitudes  are  also  of  order 
e°.  The  inconsistency  disappears  when  the  total  energy  flux  at  the  S/S 
interface  is  considered.  Recalling  that  the  energy  flux  is  proportional 
to  the  inverse  of  the  wave  speed,  the  incident  and  reflected  S  waves 
in  the  sediment  have  their  energy  flux  proportional  to  e  \  while  the 
flux  of  the  other  wave  is,  to  lowest  order,  proportional  to  e°.  The 
S  waves,  then,  carry  one  order  higher  energy  flux  than  the  other  waves 
even  though  all  the  amplitudes  are  the  same  order.  To  order  c  \  the 
incident  and  reflected  S  waves  conserve  energy  flux.  To  the  next  order, 
e°,  the  part  of  the  S  wave  energy  that  is  not  perfectly  reflected  (e^ 
corrections  to  the  amplitude)  provides  some  order  e°  energy  which  appears 
in  the  lowest  order  energy  flux  of  the  other  waves. 


We  next  consider  incident  S  and  P  waves  in  the  sediment  striking  the 
W/S  interface.  The  computed  reflection  and  transmission  coefficients  are 


p21 
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+ 

2 

2c  T  T 

21 

12  21 
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2e2T  R  „ 

21 

21  12 
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=  2e  T 
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(21) 
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for  an  incident  P  wave  and 


rP21  =~2£ 


K0>) 


COS0 


Cp21  ”  2eT21)  COS0 


rs21  =  _1 


(24) 

(25) 

(26) 


for  an  incident  S  wave.  The  lowest  order  corrections  to  Eq.  (26)  are  of 
order  e"*. 


Equations  (21)  through  (23)  give  the  same  qualitative  results  as 

Eqs.  (1)  through  (3).  Both  sets  of  equations  show  that  the  corrections 

2 

due  to  sediment  rigidity  are  of  order  e  at  the  W/S  interface  when  a  P 

wave  is  incident,  i.e.,  the  sediment  can  be  accurately  treated  as  a  fluid 

at  the  W/S  interface  whenever  a  P  wave  is  incident.  Recalling  that  the 

incident  S  wave  is  at  most  of  order  e  [Eqs.  (3),  (6),  (10)],  the  corrections 

2 

to  P  wave  amplitudes  due  to  S  wave  conversion  are  at  most  of  order  e 
[Eqs.  (24)  and  (25)].  This  leads  to  the  conclusion  that  the  sediment 
properties  at  the  W/S  interface  are  very  simple.  It  is  a  fluid  as  far 
as  P  waves  are  concerned.  The  S  wave  incident  from  below  is  perfectly 
reflected  with  a  180°  phase  shift  [Eq.  (26)]. 

The  mechanism  for  sediment  S  wave  excitation  to  lowest  order  in  e 
can  now  be  deduced  from  the  combined  action  of  the  reflection  and  trans¬ 
mission  coefficients.  As  far  as  p  waves  are  concerned,  the  sediment  is 
a  fluid  at  the  W/S  interface  and  negligible  energy  is  transferred  into  S 
waves.  An  order  e  S  wave  is  excited  in  the  sediment  by  the  P  wave 
striking  the  S/S  interface.  The  sediment  S  wave  has  no  further  effect. 

It  is  essentially  perfectly  reflected  from  the  W/S  and  S/S  interfaces 

2 

and  produces,  at  most,  order  z  corrections  to  the  P  wave  returned  to 
the  water  column.  The  P  wave  reflected  from  the  S/S  interface  does, 
however,  have  an  important  role  to  play.  It  picks  up  an  order  c 
correction  to  its  amplitude  at  the  S/S  interface.  This  correction 
travels  through  the  W/S  interface  and  affects  the  energy  returned  to 
the  water. 

The  following  picture  now  emerges.  Sediment  rigidity  affects  bottom 
reflection  loss  by  means  of  energy  transferred  to  sediment  shear  waves 
at  the  sediment-substrate  interface.  Sediment  compressional  waves  play 
a  major  role,  in  that  they  excite  the  shear  wave  at  the  sediment-substrate 
interface  and  are  responsible  for  carrying  the  influence  of  the  shear  wave 
back  through  the  sediment  and  into  the  water.  The  major  effect  of  sedi¬ 
ment  rigidity  is  then  to  provide  an  energy  sink  at  the  sediment-substrate 

interface.  We  note  that  this  situation  will  be  slightly  modified  if 
2 

higher  order  (e  and  above)  corrections  are  important. 
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The  empirical  classification  of  sediments  as  thick  or  thin  can  now 

be  understood  and  quantified.  For  a  thick  sediment  layer  the  P  wave  has 

its  turning  point  well  above  the  S/S  interface,  and  it  does  not  signifi- 

2 

cantly  interact  with  the  S/S  interface.  Since  only  order  £  corrections 
are  possible  due  to  the  refracted  P  waves  at  the  W/S  interface,  sediment 
rigidity  has  a  negligible  effect  on  RL.  In  this  case  the  sediment  can 
be  accurately  treated  as  a  fluid.  For  a  thin  sediment  layer,  the  P 
wave  interacts  significantly  with  the  S/S  interface,  where  sediment  S 
waves  are  excited.  In  this  case  the  sediment  must  be  treated  as  a  solid. 
The  question  of  whether  a  sediment  layer  is  thick  or  thin  is  then  trans¬ 
lated  into  the  more  easily  quantified  question  of  whether  the  P  wave 
interacts  with  the  sediment-substrate  interface  (thick),  or  not  (thin). 

The  reason  for  the  similarity  of  the  FF  and  SF  cases,  noted  in  the 

previous  section,  can  now  be  understood.  Equations  (21)  through  (23) 

can  also  be  applied  to  the  S/S  interface  when  the  substrate  is  taken  to 

2 

be  a  fluid;  only  order  e  corrections  are  possible.  Negligible  energy 
is  transferred  to  sediment  S  waves.  The  sediment  can  be  accurately 
treated  as  fluid  whenever  the  substrate  is  a  fluid — this  leads  to  the 
conclusion  that  a  solid  substrate  is  required  for  sediment  rigidity  to 
affect  RL.  This  statement  can  be  generalized  to  include  the  possibility 
of  S  wave  excitation  at  an  interface  within  a  layered  sediment,  provided 
the  layer  has  a  significantly  larger  shear  velocity. 

B.  Qualitative  Ray  Picture 

The  transmission  and  reflection  coefficients  derived  above  will  now 
be  applied  to  the  development  of  a  ray  model  of  the  qualitative  effects 
of  sediment  rigidity.  The  reflection  and  transmission  coefficients 
obtained  above  will  be  used  to  determine  the  ray  amplitudes  at  each 
interface.  A  single  homogeneous  sediment  layer  will  be  considered.  The 
expansion  parameter  will  be  subscripted  by  U  or  L  to  identify  the  values 
at  the  W/S  or  S/S  interfaces,  respectively.  This  allows  the  corrections, 
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due  to  sediment  rigidity,  to  be  identified  at  each  interface  and  allows 

us  to  consider  the  qualitative  effects  of  gradients  in  shear  speed  by 

setting  £>£  . 

L  u 

The  results  are  shown  schematically  in  Figs.  8,  9,  and  10.  The 
solid  lines  are  the  ray  paths  of  P  waves  and  the  dashed  lines  are  the 
ray  paths  of  S  waves.  Since  the  S  wave  speed  in  marine  sediments  is 
so  small  compared  to  the  sound  speed  in  water,  the  S  wave  rays  are 
almost  normal  to  the  plane  of  stratification.  Figure  8  shows  the 
effect  of  sediment  rigidity  if  S  wave  attenuation  is  high  enough  for  the 
S  wave  to  be  totally  absorbed  before  crossing  the  layer.  Figure  9  shows 
the  additional  ray  paths  resulting  from  the  first  reflection  of  the 
sediment  S  waves  and  Figure  10  shows  those  due  to  the  second  reflection. 

In  Figure  8  the  S  wave  attenuation  is  high  enough  for  an  S  wave 
to  be  completely  absorbed  in  traveling  one  way  through  the  sediment. 

The  incoming  water  wave  is  reflected  at  the  W/S  interface  with  corrections 

2  2 
of  order  [Eq.  (1)].  The  transmitted  P  wave  also  has  an  order 

correction  [Eq.  (2)]  which  is  neglected.  A  small  amplitude  S  wave  of 
2 

order  e  is  also  excited  [Eq.  (3)].  The  transmitted  P  wave  strikes  the 
u  2 

S/S  interface  where  an  order  e  +e  S  wave  is  excited  [Eq.  (6)].  The 

^  L  2 

reflected  P  wave  and  the  substrate  P  and  S  waves  suffer  order  EL+eL 

corrections  [Eqs.  (A),  (5),  and  (7)].  The  P  wave  carries  its  corrected 

amplitude  back  through  the  sediment  to  the  W/S  interface.  Here  a  small, 

2 

order  S  wave  is  excited  [Eq.  (23)}.  The  transmitted  P  wave  carries 

the  c  +e,  corrections  into  the  water  and  also  picks  up  a  negligible 
L*  Li 

2 

order  e  correction  from  the  W/S  interface.  The  reflected  P  wave  carries 
u  2 

corrections  of  the  order  eL+e^  back  down  through  the  sediment  where  it 
interacts  with  the  S/S  interface  and  the  process  repeats.  The  directly 
reflected  ray  and  the  rays  excited  in  the  water  by  the  sediment  P  waves 
combine  to  form  a  wave,  A  ,  which  is  of  order  t°  with  e  +e  corrections. 

O  Li  Li 

Figure  9  shows  the  additional  rays  set  up  when  the  S  wave  attenuation 

is  reduced  to  allow  the  sediment  S  waves  to  strike  one  interface  before 

2 

being  absorbed.  The  order  £u  ray  originating  at  the  W/S  interface  is 
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-  COMPRESSIONAL  WAVE 

- SHEAR  WAVE 


FIGURE  8 

THE  RAY  PICTURE:  SEDIMENT  SHEAR  WAVE 
ABSORBED  BEFORE  REFLECTION 
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—  COMPRESSIONAL  WAVE 

- CUPAP  WAVE 


FIGURE  9 

THE  RAY  PICTURE:  THE  FIRST  REFLECTION 
OF  THE  SEDIMENT  SHEAR  WAVES 
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ARL:UT 
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COMPRESSIONAL  WAVE 


FIGURE  10 

THE  RAY  PICTURE:  THE  SECOND  REFLECTION 
OF  THE  SEDIMENT  SHEAR  WAVES 


ARL:UT 
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PJV- GA 
8- 17-  79 
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2 

totally  reflected  from  the  substrate  [Eq.  (8)]  and  carries  its  e 

2  u 

magnitude.  Corrections  of  order  e  e  are  neglected.  It  sets  up  order 

2  u  L  2 
e  P  and  S  waves  in  the  substrate  and  a  reflected  P  wave  of  order  e  in 
u  2  u 

the  sediment.  The  P  wave  transmits  its  e^  magnitude  into  the  water  and 

into  a  reflected  P  wave  at  the  W/S  interface  [Eqs.  (21)  and  (22)].  The 

2 

reflected  S  wave  is  of  order  and  is  neglected  [Eq.  (23)].  The 

downward  propagating  P  wave  strikes  the  S/S  interface  where  it  sets  up 

2  2 

order  e  waves  in  the  substrate  and  an  order  e  reflected  P  wave 
u  u 

[Eqs.  (4),  (5),  and  (7)].  No  significant  S  wave  is  excited  in  the 
sediment  [Eq.  (6)].  The  upgoing  P  wave  then  continues  the  process.  The 
S  wave  originating  at  the  substrate  interface  is  totally  reflected  from 

o 

the  W/S  interface  [Eq.  (26)]  and  retains  its  order  e  +c^  magnitude.  A 
reflected  P  wave  of  order  e^e^  is  excited  [Eq.  (24)]  along  with  a 
transmitted  P  wave  of  amplitude  e^e^  in  the  water  [Eq.  (25)].  The 
reflected  P  wave  carries  its  magnitude  through  a  reflection  at  the 
S/S  interface  [Eq.  (4)]  where  it  also  excites  P  and  S  waves  of  order 
e  e  in  the  substrate  [Eqs.  (5)  and  (7)].  No  significant  sediment 

U  L| 

S  wave  is  set  up  [Eq.  (6)].  The  upward  traveling  P  wave  is  again  reflected 

from  the  W/S  interface.  The  reflected  sediment  P  wave  and  the  transmitted 

wave  in  the  water  carry  an  order  e  e.  magnitude  [Eqs.  (21)  and  (22)].  No 

u  L 

significant  S  wave  is  excited  [Eq.  (23)].  The  P  wave  then  continues  the 

process.  The  sum  of  the  rays  returned  to  the  water  sets  up  a  wave  A„ , 

2  1 

whose  magnitude  is  second  order,  i.e.,  order  e  or  e  e, . 

u  u  L 

Figure  10  shows  the  situation  for  an  attenuation  small  enough  to 

allow  the  sediment  S  waves  to  be  reflected  a  second  time  before  being 

absorbed.  The  S  wave  originating  at  the  W/S  interface  produces  negligible 

effects  when  it  strikes  the  W/S  interface.  The  reflected  S  wave  retains 
2 

its  order  e^  magnitude  and  no  significant  P  waves  are  excited  [Eqs.  (24), 
(25),  and  (26)].  The  S  wave  originating  at  the  S/S  interface  does  pro¬ 
duce  additional  corrections  to  the  energy  returned  to  the  water.  On 

2 

striking  the  S/S  interface  it  is  reflected  with  its  order  £  +e  magnitude 

Li  Li 

[Eq.  (10)].  It  excites  P  and  S  waves  in  the  substrate  and  a  reflected  P 

2 

wave  in  the  sediment,  all  of  order  e  +e  [Eqs.  (8),  (9),  and  (11)].  The 

Li  Li 

P  wave  in  the  sediment  travels  through  the  sediment  to  the  water 
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2 

interface  where  it  is  reflected,  retaining  its  order  e^+e^  [Eqs.  (21) 
and  (22)];  it  does  not  excite  a  significant  S  wave  [Eq.  (23)].  The 

downward  propagating  P  wave  is  reflected  from  the  substrate  interface, 

2  2 
retaining  its  e  +e  magnitude,  and  sets  up  order  c  +e  waves  in  the 
L  l*  L  L 

substrate  ]Eqs.  (4),  (5),  and  (7)].  The  reflected  P  wave  then  continues 

2 

the  process  and  an  S  wave  of  order  is  also  excited  in  the  sediment 
[Eq.  (6)].  This  S  wave  produces  no  further  effects  since  it  is 
absorbed  before  reaching  the  W/S  interface.  This  can  be  seen  by 
recalling  that  the  magnitude  of  this  S  wave  is  proportional  to  the 
P  wave  which  struck  the  S/S  interface.  The  magnitude  of  the  P  wave  is 
in  turn  proportional  to  the  magnitude  of  the  S  wave  which  originated  at 
the  S/S  interface  and  has  traveled  up  and  back  through  the  sediment.  The 
maximum  amplitude  of  the  newly  excited  S  wave  then  corresponds  to  the 
amplitude  of  the  twice  reflected  shear  wave.  By  hypothesis  (in  this 
example)  an  S  wave  is  absorbed  before  its  third  reflection.  Hence,  the 
new  S  wave  will  be  absorbed  before  it  reaches  the  W/S  interface.  The 

rays  returning  to  the  water  in  Fig.  10  combine  to  form  a  wave  A  ,  whose 

2  1 
magnitude  is  of  order  eL+eL-  This  wave  contributes  the  largest 

correction  to  the  A  wave  of  Fig.  8. 
o 

The  qualitative  ray  model  developed  in  Figs.  8  through  10  can  be 
used  to  understand  many  of  the  features  of  Fig.  5,  which  shows  the 
dependence  of  RL  on  frequency  for  the  36  m  turbidite  layer.  Figure  8 
applies  to  the  high  frequency  regime  where  shear  wave  attenuation  is 
large.  Figures  9  and  10  apply  at  lower  frequencies  and  explain  the 
origin  of  the  oscillatory  structure.  A  generalization  of  Figs.  8 
through  10  to  very  low  frequencies  shows  the  possibility  of  very  large 
interference  effects. 

According  to  Fig.  8  the  high  frequency  behavior  of  R  is  due  to 

energy  returning  to  the  water  in  the  Aq  wave  which  carries  the  effects 

of  S  wave  excitation  at  the  S/S  interface  through  corrections  of  order 
2 

e  and  e  .  The  magnitude  of  the  A  wave  depends  on  the  phase  difference 

L>  Li  O 

between  its  constituent  directly  reflected  component  and  the  component 
due  to  the  propagation  of  the  P  wave  through  the  sediment.  The  magnitude 
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of  R  is  the  result  of  interference  between  these  two  components  of  A  . 

o 

This  is  the  source  of  the  broad  oscillations  at  high  frequencies  seen 

in  both  the  fluid  sediment  and  solid  sediment  curves  of  Fig.  5.  The 

almost  constant  offset  of  4  dB  between  the  two  curves  is  due  to  the 
2 

order  t^+e  correction  to  the  magnitude  of  the  wave  transmitted 

into  the  water  by  the  P  wave  in  the  sediment.  The  S  wave  excited  at  the 

W/S  interface  has  a  negligible  effect  on  RL. 

Two  additional  interfering  waves  occur  in  the  water  when  the 

frequency  is  reduced  to  allow  the  S  waves  to  reflect  once  from  an 

2 

interface.  These  are  the  t  e  and  e  contributions  to  in  Fig.  9. 

u  L  u  2 

These  contributions  are  small  compared  to  the  corrections  due  to 

A  .  The  S  wave  excited  at  the  W/S  interface  results  in  a  series  of 

°  2 
rays  returning  to  the  water  with  magnitude  of  order  c^.  The  coherent 

sum  of  these  waves  has  a  phase  which  depends  on  both  P  and  S  wave 

propagation.  The  S  wave  dependence  occurs  through  a  phase  4>1=U)H/v7 

due  to  traveling  one  way  through  the  sediment.  The  phase  due  to  the 

P  wave  varies  more  slowly  with  frequency.  The  S  wave  excited  at  the  S/S 

interface  results  in  a  series  of  rays  whose  sum  is  of  order  £  e,  and 

whose  phase  also  has  the  usual  P  wave  part  and  the  S  wave  part, 

^“wH/v^.  The  S  wave  parts  of  the  phase  introduce  a  phase  shift  of  2r 

in  a  frequency  change  M This  results  in  a  rapid  oscillation  in 

RL  as  frequency  changes,  compared  to  the  P  wave  oscillations.  An 

oscillation  with  this  frequency  change  can  be  observed  near  90  Hz  in 

Fig.  5.  The  modulation  of  the  peak  heights  at  lower  frequencies  also 

has  this  recurrence  frequency. 

The  additional  set  of  rays  introduced  back  into  the  water  by  the 

second  S  wave  reflection,  the  A  wave  in  Fig.  10,  are  more  important. 

2  1 

They  have  a  magnitude  of  and  are  thus  larger  than  those  produced 

by  the  first  S  wave  reflection.  They  have  the  usual  P  wave  phase  change 
and  also  an  S  wave  related  phase  of  due  to  traveling  through 

the  sediment  and  back.  This  S  wave  phase  results  in  a  very  rapid 
oscillation  in  RL  as  the  frequency  is  changed.  The  frequency  separation 
between  maxima  is  which  is  half  of  Af^.  The  magnitude  of  A^ 
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is  larger  than  that  of  so  that  it  dominates  the  interference 
structure  of  RL  as  a  function  of  frequency.  Thus,  the  oscillation  at 
Af ^  is  stronger  in  Fig.  5  than  that  at  Af^. 

The  behavior  of  RL  at  lower  frequencies  where  S  wave  attenuation  is 

very  small  can  be  deduced  from  a  generalization  of  Figs.  8  through  10. 

The  S  waves  reflect  many  times  in  the  sediment.  Additional  waves  of 
2  2 

each  type  (e.  +t  c  e  ,e  )  are  produced  by  the  interactions  at  the 
L  L  u  ^  u  2 

W/S  and  S/S  interfaces.  Within  each  type  (e  for  example)  the  successive 
waves  have  a  relative  phase  shift  of  A<}>  due  to  the  additional  propaga¬ 
tion  distance  of  the  S  wave.  This  allows  a  powerful  interference 
phenomenon  to  develop.  The  A^  wave  will  interfere  destructively  with 
the  A^  wave  at  some  frequency  and  at  frequencies  separated  by  integer 
multiples  of  Af ^ ,  (Small  phase  shifts  due  to  effects  other  than  S  wave 

propagation  are  assumed  to  be  more  slowly  varying  with  frequency.)  At 
low  frequencies  the  magnitude  of  the  sum  of  all  the  in-phase  A^  waves 
can  result  in  an  order  e°  resultant  amplitude  and  hence  a  very  large 

cancellation  of  the  A  wave.  This  is  the  source  of  the  very  large  RL 

o 

peaks  in  Fig.  5. 

For  a  homogeneous  layer  (£  =e  =e)  Figs.  8  through  10  show  that  only 

three  kinds  of  waves  emerge  from  the  water-sediment  interface.  Figure  8 

shows  the  largest  of  these.  The  directly  reflected  ray  and  the  series 

2 

of  rays  carrying  a  magnitude  of  order  1+£+e  combine  to  form  a  single 

wave,  the  A  wave.  By  itself  the  A  wave  contains  the  basic  P  wave 
o  o 

interference  structure  and  is  responsible  for  the  gentle  oscillation 

structure  at  high  frequencies  in  Fig.  5.  The  A  wave  also  contains  a 

2  ° 

correction  to  its  magnitude  of  order  e+e  due  to  the  S  wave  excited  at 

the  substrate  interface  by  the  P  wave.  The  phase  of  the  Aq  wave  does 

not,  however,  contain  a  part  due  to  S  wave  propagation.  Thus  the  Aq 

wave  contains  an  increase  in  RL  due  to  sediment  rigidity,  but  carries 

the  basic  phase  of  the  P  wave.  The  second  kind  of  wave,  the  A.  wave, 

2  ^ 

is  formed  by  the  series  of  rays  of  order  e+e  shown  in  Fig.  10.  It  is  the 
most  important  correction  to  the  Aq  wave.  Its  magnitude  is  linear  in  e 
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to  lowest  order  whiLo  its  phase  contains  the  two-way  S  wave  phase 
The  interference  between  the  A  and  A^  waves  produces  the  rapid 
oscillations  in  RL  as  a  function  of  frequency  seen  in  Fig.  5.  The 

third  kind  of  wave,  the  Af)  wave,  is  formed  by  the  series  of  rays  of 

2  1  2 
order  c  in  Fig.  9.  Since  its  magnitude  is  of  order  c  it  produces 

only  small  corrections  compared  to  those  due  to  the  A  wave.  The 

phase  of  the  A.;  wave  contains  the  one-way  S  wave  phase,  4>^ .  The 

interference  of  this  wave  with  the  A  and  A.,  waves  produces  the 

o  2 

modulation  on  the  rapid  oscillations  in  Fig.  5. 

The  ray  model  also  predicts  that  the  A  wave  contributions 
should  die  out  more  rapidly  than  the  contributions  as  frequency 
increases.  Figures  9  and  10  show  that  the  magnitude  of  these  contri¬ 
butions  depends  on  the  attenuation  suffered  by  the  S  wave  in  traveling 
through  the  sediment.  The  A ^  wave  (Fig.  10)  depends  on  the  two-way 
total  attenuation  while  the  A.,  wave  (Fig.  9)  depends  only  on  the  one¬ 
way  total  attenuation.  As  S  wave  attenuation  increases  with  frequency, 
the  A^  wave  will  decrease  in  magnitude  more  rapidly  than  the  A?  wave. 

Finally,  the  ray  picture  shows  the  importance  of  S  wave  excitation 

by  P  wave  conversion  at  the  substrate  interface.  All  corrections  of 

order  c  occur  because  of  the  initial  interaction  of  the  P  wave  with  the 

substrate  in  Fig.  8.  The  order  e  corrections  in  the  A  wave  are  clearly 

o 

due  to  this  interaction  (Fig.  8).  The  A^  wave  is  excited  by  S  waves  set 

up  by  the  initial  P  wave  interaction  (Fig.  10).  Part  of  the  A.,  wave, 

the  c  e,  contribution  in  Fig.  9,  is  also  due  to  this  initial  P  wave.  These 

contributions  to  the  reflected  wave  would  all  decrease  in  magnitude  if  the 

P  wave  amplitude  were  to  be  diminished.  The  remaining  part  of  A  ,  that 
2 

labeled  in  Fig.  9,  while  basically  set  up  by  P  wave  conversion  at  the 
water  interface,  also  depends  on  the  P  wave  traveling  through  the  sediment 
layer.  The  same  P  wave  conditions  which  would  make  the  r  r.  wave  unimportant 

11  L 

would  also  affect  this  contribution  to  in  the  same  manner.  Thus,  it  is 
crucial  for  the  P  wave  amplitude  at  the  sediment-substrate  interface  to  be 
significant  if  sediment  rigidity  is  going  to  have  an  important  effect  on 

RL. 
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C.  Verification  of  Ray  Model 

In  this  section  the  computational  model  will  be  used  to  test  some 
of  the  qualitative  predictions  of  the  ray  model.  The  dependence  of  R 
on  frequency  produced  by  the  computational  model  will  be  used  to  extract 
the  magnitudes  of  the  A^ ,  A^ ,  and  waves.  This  is  accomplished  by 
using  the  interference  structure  to  separate  the  three  waves.  Homogeneous 
layers  will  be  treated  at  frequencies  near  15  Hz  and  a  grazing  angle  of 
20°.  Tables  II  and  III,  profiles  1  through  A,  give  the  parameters  of  the 

3 

layers  studied.  The  water  parameters  were  =  1.053  g/cm  , 

Cj  =  15A0  m/sec.  Parameter  ranges  are  chosen  to  illustrate  the  features 
of  the  ray  model  and  do  not  particularly  represent  parameter  ranges  in 
marine  sediments. 

The  reflection  coefficient  is  assumed  to  be  of  the  form 

| R |  =  A^  +  A^  cos<J>2  +  A2  cosUj+p)  ,  (27) 

where  A  ,  A, ,  and  A„  are  the  amplitudes  of  the  A  (order  1+c),  A 
o  1  2  2  o  1 

(order  e)  and  A?  (order  e)  waves,  and  p  is  a  frequency  independent 
phase  difference  between  the  A^  and  A2  waves.  This  form  for  |r|  assumes 
that  <J>^=wH/ v2  and  <{>2=<d2H/v2  contain  the  major  frequency  dependence  of  the 
phases.  This  is  true  in  our  case  since  the  frequency  range  will  be 
limited  so  that  the  phase  shifts  due  to  P  wave  propagation  and  the  inter¬ 
face  continuity  conditions  are  small  compared  to  those  due  to  S  wave 
propagation.  The  major  interference  pattern  in  |r|  will  be  due  to  a 
superposition  of  the  Aq  and  A^  waves.  The  presence  of  the  A 2  wave  will 
add  a  modulation  to  this  structure. 

A  typical  reflection  coefficient  to  be  analyzed  is  shown  in  Fig.  11. 
Four  of  the  maxima  and  minima  are  labeled  by  letters.  The  basic  oscilla¬ 
tory  structure  is  due  to  the  interference  of  the  A^  and  A^  waves.  The 
separation  between  peaks  occurs  when  <f>2  changes  by  2ir.  The  difference  in 
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ACOUSTIC  PARAMETERS  OF  THE  HOMOGENEOUS  SEDIMENT  LAYERS 


TABLE  III 


ACOUSTIC  PARAMETERS  OF  THE  SUBSTRATE  USED  WITH  THE  HOMOGENEOUS 

SEDIMENT  LAYERS 


Profile 

Number 

— 

SHIES oB 

C3 

(m/sec) 

kp^ 

(dB/m/kHz) 

V3 

(m/sec) 

k  _ 
s3 

(dB/m/kHz) 

1 

2.8 

5700 

0.152 

2900 

0.90 

2 

2.8 

5700 

0.152 

2900 

0.90 

3 

2.8 

5700 

0.152 

2900 

0.90 

4 

2.8 

5700 

0.152 

2900 

0.90 

5 

2.8 

5700 

0.035 

2900 

0.07 

b 

2.8 

1 

0.035 

2900 

0.07 

7 

2.8 

1 

0.035 

0.07 

8 

2.8 

1 

0.035 

2900 

0.07 

9 

5700 

0.035 

2900 

0.07 

10 

2.8 

5700 

0.035 

2900 

— 

11 

2.8 

5700 

— 

2900 

0.07 
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MAGNITUDE  OF  REFLECTION  COEFFICIENT 


FREQUENCY  -  Hz 


FIGURE  11 

EXAMPLE  OF  THE  DEPENDENCE  OF  THE  MAGNITUDE 
OF  THE  REFLECTION  COEFFICIENT  ON  FREQUENCY 


ARLiUT 
AS-79-1819-P 
PJV  -  GA 
8- 17.  79 


height  of  peaks  b  and  d  clearly  shows  the  effect  of  the  modulation 
produced  by  A2«  The  unknown  amplitudes  and  the  relative  phase  in 
Eq.  (27)  are  obtained  by  setting:  $1=4,2=0  at  peak  a;  <j>1=Ti/2,  <f>  =ir 
at  valley  b;  <^-11,  *f2=2n  at  peak  c;  q>1=  3tt/ 2 ,  <j>2=  3n  at  valley  d.  This 


generates  a  set  of  four  equations: 

R  =  A  +  A,  +  A-  cosif  ,  (28) 

a  o  1  2 

S  =  Ao  '  A1  ■  A2  sin^  »  (29) 

Rc  =  A0  +  Ai  "  a2  cos<p  ,  (30) 

Rd  =  Ao  “  A1  +  A2  Sin<il  *  (31) 


where  R^,  R^,  R^,  and  R^  are  the  measured  values  of  R  at  a,  b,  c,  and 
d.  These  four  equations  can  be  easily  solved  for  the  four  unknowns. 

The  first  parameter  to  be  investigated  in  the  test  of  the  ray 
model  is  the  sediment  shear  speed  v Profile  1  in  Tables  II  and  III 
was  used.  Figure  12  shows  RL  for  3  different  shear  speeds:  150  m/sec, 

70  m/sec,  and  a  fluid  case  v2  =  0  m/sec.  The  results  are  in  qualitative 
agreement  with  the  ray  model  predictions.  The  basic  level  of  RL  increases 
with  sediment  shear  speed.  The  rapid  oscillatory  structure,  due  to  the 
interference  of  the  Aq  and  A^  waves,  increases  in  magnitude  and  spreads 
out  in  frequency  as  v 2  increases.  The  modulation  due  to  the  small  A 2 
wave  is  clearly  evident  in  the  v2  =  150  m/sec  curve. 

Figure  13  shows  the  dependence  of  the  frequency  interval  between 
adjacent  maxima  of  RL,  Af,  on  sediment  shear  speed.  The  line  is  the  ray 
theory  prediction,  Af=v2/2H.  The  agreement  is  excellent. 

Figure  14  shows  the  dependence  of  the  amplitudes  of  Aq,  A^,  and 
A2  waves  on  sediment  shear  speed.  For  reference,  v2  =  150  m/sec 
corresponds  to  e=0.10.  In  agreement  with  the  ray  theory  predictions, 

A^  and  are  proportional  to  v2  for  small  shear  speeds.  The  third  set 
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FREQUENCY  -  Hi 


FIGURE  12 

REFLECTION  LOSS  versus  FREQUENCY  FOR 
A  40  m  THICK  HOMOGENEOUS  LAYER  WITH 
SEDIMENT  SHEAR  SPEED  AS  A  PARAMETER 

FOR  CURVE  (a)  Vj  0  m/sec;  (b)  Vj  -  30  m/sec; 
(c)  vj  ---  150  m/ sec 


FIGURE  14 


THE  DEPENDENCE  OF  THE  MAGNITUDES  OF  THE  A0, 
A1r  AND  A2  WAVES  ON  SEDIMENT  SHEAR  SPEED 


ARL:UT 
AS-79-1S22-P 
PJV-  OA 
8- 17-  79 


of  points  shows  the  change  of  Aq  due  to  sediment  rigidity,  i.e., 

AA  (e)hA  (0)-A  (e).  In  Fig.  1A,  A  increases  (RL  increases)  and  is 
o  o  o  o 

proportional  to  e  to  lowest  order,  in  agreement  with  the  ray  theory. 

The  lines  in  the  figure  are  drawn  to  illustrate  the  proportionality  to 
e.  Their  slope  does  not  reflect  any  theoretical  prediction.  All  three 
sets  of  points  show  the  effect  of  higher  order  corrections  at  larger  e. 

The  magnitude  of  A^  is  much  smaller  than  that  of  A^,  in  agreement  with 
the  ray  theory  predictions. 

The  dependence  of  Aq ,  A^,  and  A^  on  the  sediment  shear  wave 
attenuation  is  shown  in  Fig.  15.  Profile  2  in  Tables  II  and  III  was 
used.  The  earlier  discussion  of  the  ray  model  showed  that  the  magnitude 
of  the  A^  and  A^  waves  have  different  dependencies  on  the  S  wave  attenua¬ 
tion,  i.e.,  A^-exp(-2k^H)  and  A^-expC-k^H) ,  where  k^  is  the  imaginary 
part  of  the  S  wave  wave  number.  The  magnitude  of  the  Aq  wave  is  inde¬ 
pendent  of  S  wave  attenuation.  The  quantity  kg  in  Fig.  15  is  the  S  wave 
attenuation  in  dB/m/kHz.  This  is  related  to  k^  by  k^=kg (f /1000) /8 . 686 , 
since  we  assume  a  linear  dependence  of  k^  on  frequency.  The  values  of  kg 
greater  than  10  attenuate  the  S  wave  by  20  dB  in  a  distance  2H.  This 
effectively  eliminates  any  waves  other  than  those  shown  in  Figs.  8  through 
10  and  permits  the  ray  model  to  predict  a  clean  exponential  dependence  on 
kg  for  the  magnitudes  of  the  A^  and  A^  waves.  The  lines  in  Fig.  1  are 
"eyeball"  fits  to  the  data  obtained  from  the  computed  reflection  coeffi¬ 
cients.  The  slopes  are  0.0285  for  the  wave  and  0.0569  for  the  A^ 
wave.  They  differ  by  a  factor  of  2  as  predicted  by  the  ray  theory  and 
agree  very  well  with  the  predicted  slopes  of  0.028  and  0.056.  The 
change  in  Aq  due  to  S  waves,  AAq,  is  seen  to  be  independent  of  kg  as 
predicted. 

The  ray  model  also  predicts  the  dependence  of  A^  and  A^  on  the 

magnitude  of  the  sediment  P  wave  after  traveling  one  way  through  the 

sediment,  A  -exp(-k.H),  A. -exp(-2k.H)  and  AA  ~exp(-2k.H) .  For  P  waves 
2  11  i  o  i 

the  exponential  decay  can  come  from  two  sources:  the  attenuation  due 


A3 


MAGNITUDE 


FIGURE  15 

THE  DEPENDENCE  OF  THE  MAGNITUDES  OF  THE 
A0,  A|,  AND  A2  WAVES  ON  THE  ATTENUATION 
OF  THE  SEDIMENT  SHEAR  WAVE 


ARL:UT 
AS-79-1823-P 
PJV ■ GA 
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to  traveling  through  the  sediment  and  also  from  the  evanescent  nature 

of  the  P  wave  if  the  sediment  sound  speed  is  larger  than  that  in  water. 

Figure  16,  obtained  using  profile  3  in  Tables  II  and  III,  shows  the 

dependence  of  AA  ,  A  ,  and  A,  on  k  ,  the  P  wave  attenuation  in  dB/m/kHz. 
o  1  2  p 

The  lines  are  again  "eyeball"  fits  to  the  calculated  magnitudes.  The 

slopes  of  the  lines  are  0.059  for  A» ,  0.123  for  A,,  and  0.133  for  AA  . 

2  1  o 

The  A,  and  AA  slopes  are  nearly  double  that  of  the  A„  wave  in  agree- 
I  o  2 

meat  with  the  ray  model.  The  values  of  the  slopes  also  agree  well  with 

the  predicted  values  of  0.056  and  0.112.  The  values  of  k  were  computed 

P 

from  the  actual  imaginary  part  of  the  wave  number,  including  corrections 
due  to  grazing  angle. 


The  ray  model  strictly  predicts  no  important  S  wave  effects  for 

the  cutoff  situation,  where  the  sediment  P  wave  speed  is  larger  than 

the  sound  speed  in  water.  The  incoming  ray  in  the  water  is  specularly 

2 

reflected  with  only  order  £  corrections.  The  computational  model, 
based  on  a  wave  theory,  allows  the  penetration  of  the  P  wave  into  the 
sediment.  It  is  evanescent  and  decays  exponentially.  Generalizing  the 
ray  model  prediction  for  the  case  of  P  wave  attenuation  to  the  case  of 

the  evanescent  P  wave,  one  obtains  the  same  predictions,  with  k.  now 

2  1/2  /  1 

given  by  k^=27Tf  [  (c^cos^/c^)  -1]  /c2'  Since  this  exponential  attenua¬ 

tion  is  proportional  to  frequency,  it  can  be  written  in  terms  of  a 
frequency  independent  decay  constant  k^,  in  units  of  dB/m/kHz,  i.e., 

k  =k'f/8686. 
i  P 


Figure  17  shows  the  dependence  of  AA  ,  A,,  and  A„  on  k'.  The 

o  1  2  p 

actual  parameter  varied  in  this  case  was  the  P  wave  velocity  in  the 

sediment.  Profile  4  in  Tables  II  and  III  was  used.  The  lines  are  "eyeball" 

fits  to  the  calculated  magnitudes.  The  measured  slopes  are  0.162  for  AAq, 

0.162  for  A^,  and  0.072  for  A2;  the  slope  of  the  A^  line  is  approximately 

half  that  of  the  slopes  of  the  AA  and  A  lines  in  agreement  with  the  ray 

o  1 

model  predictions.  The  magnitude  of  the  slopes  also  agrees  well  with  the 
predicted  values  of  0.070  and  0.140. 
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MAGNITUDE 


COMPRESSIONS.  WAVE  ATTENUATION  -  dBWkHi 


FIGURE  16 

THE  DEPENDENCE  OF  THE  MAGNITUDES  OF  THE 
A0,  A,,  AND  A2  WAVES  ON  THE  ATTENUATION 
OF  THE  SEDIMENT  COMPRESSIONAL  WAVE 
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MAGNITUDE 


One  ot  tin*  Important  com- Ins  ions  of  tills  paper  is  that  sediment 

rigidity  can  be  neglected  at  the  water-sediment  interlace.  The  ray 

model  supports  tills  conclusion,  since  it  predicts  that  the  lowest  order 

effects  depending  on  the  S  wave  velocity  at  the  W/S  interface  are 
2 

contained  in  the  >  and  r  i  .  waves  seen  in  Fig.  9.  Since  tsO.l.  these 
u  u  1.  ’ 

are  small  corrections  compared  to  the  i  corrections  seen  in  Figs.  8 
and  10  due  to  P  wave  conversion  at  the  S/S  interface. 

Figure  18  illustrates  the  relative  Insensitivity  of  R  to  the 

sediment  S  wave  velocity  at  the  W/S  Interface  (v,  )  compared  to  that  at 

2u 

the  S/S  interface  (v.;]).  Figures  18(a)  and  18(h)  use  the  computational 

model  to  tollow  the  development  of  | R| ,  starting  with  a  homogeneous 

layer  ill  which  v.,  *=  100  m/sec  and  ending  with  v.,  =  100  m/sec.  To 

eliminate  interference  effects,  a  frequency  in  the  high  frequency  regime 

was  chosen.  In  Fig.  18(a)  v.,^  is  constant  while  v.  Increases.  A 

constant  gradient  is  used  which  has  the  value  ol  I /sec  when  v.,j  ~  100  m/sec. 

Three  grazing  angles  are  shown.  |r|  decreases  almost  linearly  with 

v  at  each  grazing  angle.  The  decrease  is  significant,  '102.  at  0=10°  and 

SO",  at  0*,I0".  Figure  18(h)  shows  the  additional  change  in  |l\|  which 

occurs  when  v  Is  held  constant  at  TOO  m/sec  and  v,  increased  to  lorn) 

21.  2u 

a  homogeneous  layer  with  v,  ™  100  m/sec.  As  expected,  the  factor  of 
three  change  in  v.^  now  produces  only  small  changes  in  | R | .  Compared  to 
| R |  for  the  original  100  m /sec  homogeneous  layer,  the  changes  arc  1.6% 
at  i 10°  and  /»%  at  i)"10o.  These  are  small  changes  compared  to  those 
produced  by  changing  v.,j  .  This  result  reinforces  the  conclusion  that 
Rl.  Is  relatively  Insensitive  to  the  surficlal  S  wave  velocity.  In 
fact,  fairly  accurate  values  of  [ R [  can  be  obtained  in  the  high  frequency 
regime  by  assuming  a  homogeneous  layer  with  v.,  given  by  the  S  wave 
velocity  at  the  sediment-substrate  interlace.  For  the  same  gradients  used 
in  Fig.  18(a)  the  error  produced  by  this  assumption  is  less  than  12%  over 
the  entire  range  of  gradients  and  for  grazing  angles  less  than  SO".  This 
leads  to  the  conclusion  that  v  is  the  important  S  wave  parameter  in 
the  high  frequency  regime;  gradients  and  v.,  are  not  major  factors 
Inf  luencing  Rl.. 
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MAGNITUDE  OF  REFLECTION  COEFFICIENT 


(a)  (b) 


FIGURE  18 

THE  DEPENDENCE  OF  THE  MAGNITUDE  OF  THE  REFLECTION 
COEFFICIENT  ON  THE  SEDIMENT  SHEAR  WAVE  SPEED  AT 
INTERFACES  FOR  GRAZING  ANGLES  OF  10°  (CIRCLES), 

20°  (TRIANGLES),  AND  30°  (SQUARES) 
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IV.  PARAMETER  STUDIES 


In  this  section  the  effect  of  bottom  parameter  variations  on  the 
bottom  reflection  loss  of  a  solid  sediment  will  be  studied.  The  com¬ 
putational  model  will  be  used  to  obtain  R.  The  study  will  be 
restricted  to  thin  sediment  layers  for  which  sediment  rigidity  is 
important.  Current  geophysical  data  will  be  used  to  limit  the  parameter 
ranges.  We  first  consider  homogeneous  layers  and  use  the  amplitudes  of 

the  interfering  waves  in  the  ray  picture  (A  ,  A,,  A„)  to  quantify  the 

o  1  2 

changes  in  RL  due  to  sediment  S  wave  propagation.  Next  we  treat  layers 
with  gradients,  using  mainly  the  change  in  RL  to  identify  significant 
parameter^.  Finally,  the  question  of  the  accuracy  required  in  sediment 
parameter  values  to  accurately  predict  RL  is  addressed.  This  is  done 
by  varying  the  parameters  of  a  typical  deep  sea  sediment. 

A.  Homogeneous  Layers 

The  ray  picture  presented  earlier  shows  that  RL  is  due  to  the 

interference  between  the  A  ,  A, ,  and  A.  waves.  The  interference 

o  1  L 

structure  can  be  quite  strong  at  low  frequencies  where  sediment  S  wave 
attenuation  is  small.  At  high  frequencies  the  increased  S  wave  attenu¬ 
ation  effectively  eliminates  the  A^  and  A^  waves  and  destroys  the 

interference  pattern.  The  A  wave  remains  to  carry  the  effects  of 

o 

sediment  rigidity  at  high  frequencies.  The  A^  wave  also  carries  the 

effects  due  to  a  fluid  sediment.  The  difference  AA  between  the  A 

o  o 

wave  for  fluid  and  solid  sediments  is  entirely  due  to  sediment  rigidity. 

In  this  section  the  results  of  our  study  of  the  influence  of 
parameter  variations  in  a  homogeneous  sediment  layer  on  AAq,  A^,  and 
A^  will  be  presented.  In  these  studies  the  grazing  angle  was  held  at 
20°,  and  the  wave  magnitudes  were  obtained  at  a  frequency  of  14  Hz. 

The  changes  produced  by  variations  of  a  single  parameter,  over  the 
range  given  by  geophysical  data,  will  be  presented.  The  bottom 
structure  is  given  in  Tables  II  and  I’ I:  Table  II  gives  sediment 
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parameters,  and  Table  III  gives  substrate  parameters.  In  all  cases 

the  sound  speed  in  water  was  taken  to  be  1540  m/sec  and  the  water 

3 

density  to  be  1.053  g/cm  .  First,  the  effect,  of  parameter  variation 
on  AAq  will  be  discussed;  its  behavior  will  hold  in  both  high  and  low 
frequency  regimes.  Next,  we  will  discuss  the  behavior  of  the  A  and 
A2  waves  which  are  important  in  the  low  frequency  regime.  Three 
categories  of  parameters  will  emerge:  those  that  are  definitely 
important,  those  which  can  be  neglected,  and  those  which  are  important 
under  certain  conditions.  The  criterion  used  is  that  a  parameter  is 
important  if  it  produces  a  change  of  10%  when  varied  over  half  the 
possible  range  of  variation  given  by  geophysical  data.  The  parameter 
ranges  are  given  in  Table  IV.  We  restrict  ourselves  to  high  porosity, 
deep  sea  sediments  for  which  our  computational  model  applies.  The 
profiles,  from  Tables  II  and  III,  used  in  each  study  are  also  given  in 
Table  IV. 


As  expected  from  the  ray  picture  and  the  expansion  theory  of  the 

reflection  and  transmission  coefficients,  the  sediment  shear  speed 

greatly  affects  AAq.  Not  only  is  a  finite  shear  speed  necessary  to 

describe  thin  sediments,  but  changing  it  by  50  m/sec  produces  a  50% 

increase  in  AAq.  Surprisingly,  the  only  other  important  parameters 

are  the  sediment  P  wave  speed  and  density.  This  is  surprising  since 

for  the  corresponding  fluid  sediments  we  find  that  the  variation  of 

these  parameters  produces  very  small  changes  in  RL  (<4%) .  For  the  solid 

sediment  case,  a  change  in  P  wave  speed  of  80  m/sec  produced  50 '1  changes 

3- 

in  AA  and  a  change  of  0.25  g/cm  in  sediment  density  changed  AA  by 
about  15%. 


In  the  category  of  negligible  parameters  for  AAq  are  substrate 
density  and  substrate  P  and  S  wave  attenuations.  These  parameters 
only  enter  through  the  boundary  conditions  at  the  sediment-substrate 
interface.  The  negligible  effect  of  the  attenuations  is  reasi table 
since  they  are  only  a  small  part  of  the  total  complex  wave  velocities. 
Sediment  P  wave  attenuation  can  also  be  ignored.  This  is  due  to  the 
low  frequency,  and  to  sediment  thickness. 
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TABLE  IV 


ACOUSTIC  PARAMETER  RANGES,  PROFILE  NUMBERS  IN  TABLES  II  AND  III 
USED  IN  PARAMETER  STUDIES,  AND  REFERENCE  NUMBERS 


Parameter 

Range 

Profile 

Number 

Reference 

V2 

m/sec 

50-150 

1 

10 

ks2 

dB/m/kHz 

8.7-26.1 

2 

11 

k2 

dB/m/kHz 

0-0.2 

3 

18 

C2 

m/ sec 

1500-1620 

mm 

19 

P2 

g/cm3 

1.1-1. 6 

6 

20 

V3 

m/sec 

1700-3600 

7 

21 

C3 

m/sec 

3500-6700 

8 

21 

P3 

g/cm3 

2. 1-3.0 

9 

21 

k  o 

s3 

dB/m/kHz 

0.06-0.50 

10 

11 

k  . 

P3 

dB/m/kHz 

0.02-0.05 

11 

18 
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The  category  of  conditionally  important  parameters  for  AA 

o 

contains  the  sediment  S  wave  attenuation  and  the  substrate  P  and  S  wave 
velocities.  As  expected  from  Fig.  15,  varying  the  shear  wave  attenua¬ 
tion  over  most  of  its  range  has  an  insignificant  effect  on  AA  .  Only 

o 

for  very  low  values  of  k  .  is  AA  affected.  This  is  probably  due  to 

s2  o  ^ 

the  extreme  interference  effects  for  small  k  _  producing  many  order  eA 

s2 

rays,  which  then  combine  into  an  order  e  contribution.  These  additional 

contributions  would  not  have  the  expected  e  scaling  and  would  appear  as 

part  of  AA^  in  our  analysis  of  j  R | .  The  effect  of  substrate  compressional 

speed  on  AAq  is  shown  in  Fig.  19.  For  large  c^,  between  5000  and 

6600  m/sec,  the  change  in  AA  is  about  10%.  For  small  c  .  between 

o  3 

3400  and  5000  m/sec,  the  change  in  AA^  is  significant,  particularly 
at  the  lower  values.  The  substrate  P  wave  speed  is  then  a  significant 
parameter  for  small  values,  but  is  negligible  if  its  value  is  large 
enough.  The  same  is  true  of  the  substrate  S  wave  speed,  the  separation 
between  the  two  regimes  being  at  2600  m/sec. 

One  additional  parameter  is  the  sediment  thickness.  For  homogeneous 

layers  the  differential  equations  for  sediment  P  and  S  waves  can  be 

written  in  terms  of  scaled  depth  Z=z/H,  where  z  is  the  depth  below  the 

water-sediment  interface.  The  complex  wave  velocities  are  then 

scaled  to  the  velocity  V=wH.  The  interface  continuity  conditions  can 

also  be  written  in  terms  of  the  same  velocity,  V.  This  means  that  for 

homogeneous  layers  the  frequency  and  depth  en^er  only  through  the 

parameter  V.  The  dependence  of  RL  on  depth  is  then  the  same  as  the 

dependence  on  frequency.  In  fact,  RL  obtained  for  u)*au  and  H=H  is 

o  o 

the  same  as  that  for  w=w0  and  H^aH^,  i.e.,  increasing  the  frequency 
is  the  same  as  increasing  the  layer  thickness.  This  relation  between 
id  and  H  will  hold  only  approximately  for  layers  with  gradients.  In 
this  case  there  are  additional  scale  lengths  in  the  problem,  those 
associated  with  the  gradients  which  prevent  the  differential  equations 
from  being  scaled. 

The  dependence  of  A^  and  A^  on  the  bottom  parameters  has  a  slightly 
different  structure.  The  sediment  density  shifts  from  being  an  important 
parameter  to  being  a  negligible  parameter  whereas  the  sediment  S  wave 
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attenuation  shifts  from  the  conditional  to  important  category.  All  other 
parameters  stay  in  the  same  categories. 

B.  Layers  with  Gradients 

This  section  presents  the  effect  of  parameter  gradients  on  RL  of  a 

solid,  deep  sea  sediment.  This  study  of  gradient  effects  is  accomplished 

by  allowing  one  parameter  at  a  time  to  have  a  gradient  in  an  otherwise 

homogeneous  layer.  The  parameters  of  the  homogeneous  layer  are  typical 

of  deep  sea  sediments.  Using  the  subscript  convention  of  the  expansion 

3 

theory  above,  the  parameters  are  P.  “  1.053  g/cm  ,  c  =  1540  m/sec, 

3 

P2  *  1.25  g/cm  ,  c2  -  1525  m/sec,  k^2  =  0.100  dB/m/kHz,  v?  =  100  m/sec, 
kg2  =  10.0  dB/m/kHz,  p^  =  2.8  g/cm^',  c^  =  5700  m/sec,  k  ^  =  0.035  dB/m/kHz, 
v  =  2900  m/s ec,  and  k  =  0.070  dB/m/kHz.  The  layer  thickness  is  H  =  40  m. 

j  S  j 

The  effect  of  gradients  on  the  high  frequency  regime  was  obtained  by 
examining  the  change  in  RL  as  a  function  of  grazing  angle  with  the  gradient 
as  a  parameter  at  a  fixed  frequency  (50  Hz).  The  effect  on  the  low 
frequency  interference  regime  was  obtained  in  a  similar  manner,  with  the 
frequency  chosen  to  be  one  of  the  peaks  in  RL.  For  most  parameters  this 
frequency  is  19.45  Hz.  For  the  case  of  the  gradient  in  sediment  S  wave 
velocity  the  frequency  was  chosen  for  each  gradient  to  be  the  largest 
interference  peak  near  the  original  frequency  of  19.45  Hz.  This 
frequency  depends  on  the  gradient  since  the  shear  wave  phase  accumulated 
in  traveling  through  the  sediment  depends  on  the  actual  gradient. 

The  first  parameter  to  be  discussed  is  the  S  wave  speed.  Typical 
gradients  of  v 2  are  expected  to  be  about  4,65  sec  ^  (Ref.  10). 

Figures  20  and  21  show  the  effect  of  the  gradient  of  S  wave  speed  on 
RL.  At  50  Hz  (Fig.  20),  the  high  frequency  regime,  RL  increases 
steadily  with  increasing  gradient.  This  is  the  result  to  be  expected 
from  the  ray  picture.  The  major  effect  of  the  gradient  at  high  fre¬ 
quencies  is  to  increase  the  S  wave  speed  at  the  lower  interface,  thus 
increasing  the  sediment  S  wave  amplitude.  The  increase  is  quite  large 
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with  a  peak  of  6.2  dB  at  10°,  and  decreases  at  both  higher  and  lower 
grazing  angles.  At  low  frequencies,  where  S  wave  propagation  is 
important,  the  structure  of  large  peaks  (Fig.  5)  changes  with  the 
gradient  of  S  wave  speed.  The  separation  between  peaks  increases 
from  0.65  Hz  for  the  homogeneous  layer  to  1.10  Hz  for  a  gradient  of 
5  sec  In  addition,  Fig.  21  shows  that  the  magnitude  of  the 
interference  effect  also  increases  for  larger  gradients.  The  magni¬ 
tude  of  the  increase  in  RL  is  only  about  1  dB  larger  than  that  for  the 
high  frequency  regime  of  Fig.  20.  The  increase  with  gradient  is  not  as 
regular  as  that  in  Fig.  20.  These  changes  can  be  qualitatively 
explained  by  the  ray  picture.  The  basic  effect  of  the  gradient  on  the 
change  in  RL  occurs  through  the  increase  of  S  wave  speed  at  the 
sediment-substrate  interface.  The  magnitude  of  AA  will  increase  in 
the  same  manner  in  both  high  and  low  frequency  regimes,  and  thus  the 
major  part  of  the  increase  in  RL  in  both  frequency  regimes  is 
accounted  for.  The  more  erratic  increases  at  low  frequencies  (Fig.  21) 
are  probably  due  to  the  combination  of  higher  order  corrections  to  the 
reflection  and  transmission  coefficients  (e  =  0.2  at  a  gradient  of  5) 

L 

and  the  phase  shift  between  the  A^  and  A^  waves. 

The  effect  of  the  depth  dependence  of  the  S  wave  attenuation 

is  shown  in  Fig.  22.  The  frequency  is  19. 45  Hz,  which  corresponds  to 

a  maximum  in  RL  at  0=20°.  The  depth  dependence  was  calculated  by 

assuming  that  the  ratio  between  the  logarithmic  decrements  (ratio  of 

the  imaginary  part  of  the  wave  numbers  to  the  real  part)  of  shear  and 

compressional  waves  is  constant  with  depth.  At  the  sediment  surface 

this  ratio  is  typically 

to  attenuation  gives  k^ 

depth  dependence  occurs  because  of  the  increase  of  S  wave  speed  v ^  with 

depth,  resulting  in  a  decrease  in  attenuation  with  depth.  Approximating 

this  depth  dependence  by  a  constant  gradient  and  using  parameters  typical 

of  deep  sea  sediments^’^  gives  dk  „/d  =  -0.15  dB/m  /kHz.  Figure  22 

sz/  z 


A  /A  =0.3  for  deep  sea  sediments, 
s  p 


11 


Converting 

(z)=k  2(z)[c2(z)/v2(z)]As2/A  .  The  major 
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FIGURE  22 

THE  CHANGE  IN  REFLECTION  LOSS 
versus  GRAZING  ANGLE  AT  19.45  Hz 

THE  GRADIENT  OF  SHEAR  WAVE  ATTENUATION 
(IN  UNITS  OF  dB/m2/kHi)  AS  PARAMETER 


shows  that  this  decrease  in  S  wave  attenuation  produces  a  significant 
change  in  RL  in  the  low  frequency  regime.  This  is  expected  since  a 
decrease  in  attenuation  will  increase  the  effectiveness  of  interference 
effects.  In  the  high  frequency  regime  the  gradient  in  S  wave  attenua¬ 
tion  produces  insignificant  (0.15  dB)  changes  in  RL.  This  also  is 
expected  since  in  the  high  frequency  regime  the  S  waves  are  completely 
absorbed  within  the  sediment.  The  change  in  total  attenuation  pro¬ 
duced  by  the  attenuation  gradient  is  too  small  to  keep  the  S  wave  from 
being  absorbed.  The  effect  of  the  gradient  could  at  most  increase  the 
frequency  characterizing  the  separation  of  the  high  and  low  frequency 
regimes. 

20  3 

Typical  gradients  in  density  (<0.000025  g/cm  /m)  and  P  wave 
18  ~U  2 

attenuation  (-2.5x10  dB/kH2/m  )  produced  negligible  (<0.05  dB) 
changes  in  RL  at  low  grazing  angles  (0  to  50°). 

To  study  the  effect  of  the  gradient  in  P  wave  speed  it  was 

necessary  to  remove  the  excitation  of  Stoneley  waves  at  the  sediment- 

substrate  interface.  This  was  done  by  changing  the  substrate  parameters 
3 

to  =  2.6  g/cm  ,  c^  =  4900  m/sec,  v ^  =  2400  m/s ec.  The  changes  in  RL 

shown  in  Figs.  23  and  24  are  then  due  to  the  gradient  alone  and  are  not 

related  to  the  excitation  of  interface  waves.  Gradients  of  P  wave 

22  -1 

velocity  are  typically  less  than  2  sec 

Figures  23  and  24  show  that  RL  is  a  solid  sediment  surprisingly 
sensitive  to  the  gradient  of  P  wave  speed.  At  low  frequencies  (Fig.  23) 
an  increase  of  over  4  dB  was  observed  at  low  grazing  angles  for  a 

gradient  of  2  sec  The  increase  in  RL  is  even  larger  in  the  high 

frequency  regime.  Figure  24  shows  an  increase  of  nearly  10  dB.  The 
gradient  in  P  wave  speed  is  clearly  an  important  sediment  parameter. 

C.  Sensitivity  of  RL  to  Parameter  Accuracy 

In  this  section  we  investigate  the  question  of  how  a  small  change 
in  an  ocean  bottom  parameter  influences  the  bottom  reflection  loss.  This 
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FIGURE  23 

THE  CHANGE  IN  REFLECTION  LOSS 
versus  GRAZING  ANGLE  AT  19.45  Hz 

THE  GRADIENT  OF  THE  COMPRESSIONAL  WAVE 
SPEED  (UNITS  OF  sec1)  AS  PARAMETER 
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FIGURE  24 

THE  CHANGE  IN  REFLECTION  LOSS 
versus  GRAZING  ANGLE  AT  50  Hi 
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is  done  by  perturbing  the  parameters  of  a  thin  sediment  layer  having 
realistic  parameters  and  computing  the  resulting  change  in  RL.  Parameters 
will  be  varied  one  at  a  time  so  that  their  effects  can  be  isolated.  Param¬ 
eters  will  be  grouped  into  three  categories:  high  sensitivity,  normal 
sensitivity,  and  low  sensitivity.  For  high  sensitivity  parameters  the 
fractional  change  in  RL  will  be  greater  than  the  fractional  change  in  the 
parameter.  Normal  sensitivity  parameters  will  induce  changes  in  RL  of 
the  same  size  or  slightly  less  than  the  fractional  change  in  the  parameter. 
Low  sensitivity  parameters  will  induce  a  change  in  RL  an  order  of  magnitude 
or  more  smaller  than  the  change  in  the  parameter.  Low  grazing  angles  and 
the  high  frequency  regime  will  be  treated. 

The  parameters  of  the  layer  studied  are  given  in  Table  V.  They  were 
obtained  from  the  geophysical  literature  and  are  typical  of  high  porosity, 
deep  sea  sediment  types.  Figure  25  shows  the  RL  for  0  to  45°  grazing 
angles  in  the  high  frequency  regime  at  50  Hz.  This  figure  provides 
a  reference  level  to  put  the  change  in  RL  induced  by  parameter  variations 
into  perspective.  The  expected  increase  in  RL  due  to  sediment  rigidity 
(-5  dB)  is  seen.  Bottom  parameters  will  be  changed  by  +  10%  about  the 
nominal  values  of  Table  V.  Parameters  were  put  into  the  high  sensitivity 
category  if  they  induced  a  change  in  RL  of  more  than  20%  (about  2  dB) . 

If  the  change  was  less  than  5%  (about  0.4  dB) ,  they  were  put  into  the  low 
sensitivity  category.  The  normal  category  contained  parameters  inducing 
changes  in  RL  of  between  5%  and  20%,  about  the  same  size  as  the  parameter 
change. 

The  low  sensitivity  category  contained  all  attenuations  and  their 
gradients  as  well  as  the  gradients  of  sediment  density  and  P  wave 
velocity.  The  sediment  surficial  S  wave  speed  is,  surprisingly,  also  a 
low  sensitivity  parameter.  All  were  varied  by  +10%.  In  order  of  decreas¬ 
ing  sensitivity  they  were:  V2(0.30  dB),  1^(0. 15  dB),  ^(0.08  dB), 

k  ,(0.06  dB),  p ’ (0.03  dB) ,  k'(0.02  dB) ,  and  k  ,,  k  ,,  k’  (all  inducing 
SZ  Z  s  Z  sj  pj  pZ 

changes  of  less  than  0.01  dB) .  The  prime  denoted  the  derivative  with 
respect  to  depth. 


i  i 
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FIGURE  25 

REFLECTION  LOSS  versus  GRAZING  ANGLE 
FOR  A  TYPICAL  DEEP-SEA  SEDIMENT  LAYER 
40  m  THICK  AT  50  Hi 
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TABLE  V 


ACOUSTIC  PARAMETERS  OF  THE  OCEAN  BOTTOM  USED  TO  STUDY  THE  SENSITIVITY  OF 
REFLECTION  LOSS  TO  PARAMETER  ACCURACY.  CONSTANT  GRADIENTS  ARE  ASSUMED. 


Depth 

(m) 

c 

(m/sec) 

b 

P 

(dB/m/kHz) 

P  3 

(g/cin  ) 

V 

(m/sec) 

k 

s 

(dB/m/kHz) 

water 

1530 

— 

1.03 

— 

— 

0 

■ 

116 

15.00 

40 

1 

302 

9.00 

Substrate 

4930 

0.030 

2.600 

2460 

0.030 
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Normal  sensitivity  parameters  in  order  of  decreasing  sensitivity 
were  found  to  be:  H(1.7  dB) ,  P 2 (1 . 5  dB),  P3<0.76  dB),  ^(0.68  dB), 
P^CO.13  dB  for  2%  change),  and  v^CO^S  dB) .  The  sediment  S  wave 
velocity  gradient  appears  in  this  group.  The  rather  large  changes 
induced  by  and  c^  indicate  the  importance  of  the  substrate  interface. 
The  effect  of  changing  the  sediment  thickness  by  +10%  is  seen  in 
Fig.  26.  The  large  peaks  between  20°  and  30°  are  the  result  of  changes 
in  the  level  of  RL  and  are  not  due  to  a  shift  in  the  RL  curve  with 
angle.  Figure  27  shows  the  effect  of  changing  the  sediment  surficial 
density  by  ±10%.  The  induced  changes  are  also  concentrated  at  low 
grazing  angles. 

The  high  sensitivity  parameters  are  c, ,  c^,  and  v^.  Varying  the 
sound  speed  in  water  by  +1%  induced  the  change  in  RL  shown  in  Fig.  28. 

The  change  in  RL  for  -1%  is  nearly  the  mirror  image  of  that  for  +1%. 

The  changes  are  most  important  at  very  low  grazing  angles  (~5°)  and  are 
small  above  15°.  Inspection  of  the  RL  versus  0  curve  shows  that  the 
changes  near  5°  are  not  due  to  a  slight  shift  in  angle  of  the  rapidly 
increasing  RL  at  low  angles,  but  rather  are  due  to  an  actual  change  in 
RL.  The  change  in  RL  induced  by  changing  by  +10%  is  shown  in  Fig.  29. 
The  corresponding  changes  for  a  fluid  sediment  are  also  shown.  The 
changes  are  not  symmetric  for  +10%  changes  in  c^.  For  the  solid 
sediment  RL  decreases  quite  substantially  at  small  angles.  For  the 
fluid  sediment  case  the  changes  are  much  smaller  except  that  the  large 
peak  occurs  near  20°  in  the  +10%  case;  this  peak  is  due  to  the  excitation 
of  a  Stoneley  wave  at  the  substrate  interface  by  the  evanescent  P  wave 
in  the  sediment.  Note  that  this  does  not  occur  in  the  solid  sediment 
case,  i.e.,  sediment  rigidity  greatly  reduces  the  coupling  into  interface 
waves  at  the  substrate  interface.  The  very  large  sensitivity  to  c^  is 
reduced  substantially  if  the  ratio  c2^c\  ^"s  constant*  Varying  c^ 

(and  hence  c^)  by  1%  then  produces  changes  of  at  most  0.28  dB,  an  order 
of  magnitude  smaller  than  the  2.7  dB  change  seen  in  Fig.  28.  The  effect 
of  varying  substrate  S  wave  velocity  is  given  in  Fig.  30.  The  changes 
are  nearly  symmetric  for  ±10%  changes.  The  large  increase  near  45°  is 
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CHANGE  IN  REFLECTION  LOSS  -  dB 


CHANGE  IN  REFLECTION  LOSS  -  dB 


FIGURE  29 

THE  CHANGE  IN  REFLECTION  LOSS  versus 
GRAZING  ANGLE  FOR  ±10%  CHANGES 
IN  THE  SURFICIAL  COMPRESSIONAL 
WAVE  SPEED  IN  THE  SEDIMENT 


ARL:UT 
AS-79-1837-P 
PJV- GA 
8-17-  79 
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CHANGE  IN  REFLECTION  LOSS  -  dB 


due  to  a  slight  shift,  to  lower  grazing  angles,  of  the  rapidly  increasing 
RL  near  50°  due  to  substrate  shear  wave  excitation  near  its  critical 
angle. 


D.  Additional  Topics 

Two  final  topics  are  considered  in  this  section.  First  is  the 
effect  of  sediment  rigidity  on  the  excitation  of  interface  waves  on  the 
sediment-substrate  interface.  For  the  case  of  a  fluid  sediment,  Stoneley 
waves  traveling  along  the  substrate  interface  have  been  identified  as  a 
possible  loss  mechanism  at  low  grazing  angles,  producing  very  large  RL 
over  a  very  small  interval  of  grazing  angles.  The  characteristics  of 
these  interface  waves  for  fluid  sediments  have  been  studied  in  detail. 

The  last  topic  is  a  discussion  of  sediment  rigidity  effects  in  sands,  an 
important  sediment  type  in  shallow  water,  continental  shelf  areas. 

The  bottom  parameters  used  in  the  study  of  sediment  rigidity  on 

RL  were  taken  from  the  work  of  Hawker  so  that  the  Stoneley  wave  peak 

in  RL  for  the  fluid  sediment  case  could  be  unambiguously  r'dentif ied .  . 

3 

The  parameters  of  the  water  column  are:  =  1.053  g/cm  ,  c^  =  1540  m/sec. 

The  sediment  layer  is  100  m  thick.  Sediment  parameters  at  the 

3 

water-sediment  interface  are:  =  1.27  g/cm  ,  c^  =  1527  m/sec, 

k  -  =  0.057  dB/m/kHz,  k  „  =  8.0  dB/m/kHz.  The  sediment  shear  wave 
p2  s2 

speed  was  varied  to  determine  the  effect  of  sediment  rigidity.  A 

homogeneous  layer  was  assumed  except  for  a  constant  sound  speed  gradient, 

-1  3 

=  1.2  sec  .  The  substrate  parameters  are:  =  2.6  g/cm  , 

c„  =  5700  m/sec,  k  „  =  0.03  dB/m/kHz,  v_  =  2900  m/sec,  k  =  0.030  dB/m/kHz. 
3  p3  3  s3 

Figures  31  through  33  show  the  effect  of  sediment  rigidity  on 
interface  wave  excitation  at  frequencies  of  25,  50,  and  100  Hz, 
respectively.  The  dashed  curves  were  obtained  for  the  fluid  sediment 
case.  The  peaks  in  RL  near  18°  are  due  to  Stoneley  wave  excitation. 

The  solid  curves  in  each  figure  show  the  effect  of  sediment  rigidity  for 
”  100,  200,  and  400  m/sec.  A  value  of  400  m/sec  is  typical  at  a 
depth  of  100  m  for  a  clay  sediment. 


REFLECTION  LOSS  -  dB 


GRAZING  ANGLE  -  deg 


FIGURE  31 

REFLECTION  LOSS  versus  GRAZING  ANGLE 
FOR  A  100  m  THICK  CLAY  LAYER  AT  25  Hr 

THE  PARAMETER  IS  THE  SEDIMENT  SHEAR 
WAVE  SPEED  IN  UNITS  OF  m/*ec 


ARL  :UT 
AS-79-1839-P 
PJV- GA 
8-17- 79 
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GRAZING  ANGLE  -  deg 


FIGURE  32 

REFLECTION  LOSS  versus  GRAZING  ANGLE 
FOR  A  100  m  THICK  CLAY  LAYER  AT  50  Hz 

THE  PARAMETER  IS  THE  SEDIMENT  SHEAR 
WAVE  SPEED  IN  UNITS  OF  m/sec 


ARL:UT 
AS-79-1840-P 
PJV-  GA 
8-  17-  79 


GRAZING  ANGLE  -  deg 


FIGURE  33 

REFLECTION  LOSS  versus  GRAZING  ANGLE 
FOR  A  100  m  THICK  CLAY  LAYER  AT  100  Hr 

THE  PARAMETER  IS  THE  SEDIMENT  SHEAR 
WAVE  SPEED  IN  UNITS  OF  m/sec 


ARL.UT 
AS-79-1841-P 
PJV-  GA 
8-  17-  79 


At  25  Hz  (Fig.  31)  and  50  Hz  (Fig.  32),  a  small  increase  in  the 
shear  speed  (v^  =  100  m/sec)  greatly  increases  the  energy  loss  to  the 
interface  wave.  Increasing  the  S  wave  velocity  further  shifts  the 
peak  location  to  higher  angles.  The  peak  is  both  broadened  and  reduced 
in  magnitude.  At  v^  =  400  m/sec  the  peak  at  25  Hz  can  still  be  dis¬ 
tinguished  although  it  is  now  very  broad.  At  50  Hz  the  peak  can  no 
longer  be  unambiguously  identified  but  appears  as  part  of  the  set  of 
P  wave  interference  peaks.  A  general  increase  in  RL  above  about  17° 
is  seen  for  all  values  of  S  wave  speed.  This  is  expected  since  the 
location  of  the  Stoneley  wave  (v^O)  in  angle  occurs  very  near  the 
critical  angle  below  which  the  P  wave  interaction  with  the  substrate  is 
negligible.  For  larger  angles  sediment  S  waves  are  excited  at  the 
substrate  interface  and  RL  is  increased. 

At  100  Hz  (Fig.  33)  the  energy  loss  to  interface  waves  does  not 
greatly  increase  for  v^  =  100  m/sec.  Rather,  the  interface  wave  peak 
is  broadened  and  lowered  in  magnitude.  For  v0  =  200  m/sec  and  400  m/sec 
the  existence  of  an  interface  wave  peak  in  RL  cannot  be  clearly  observed. 
The  RL  shows  the  expected  increase  due  to  sediment  S  wave  excitation  in 
angles  above  about  18°. 

One  possible  explanation  for  the  results  in  Figs.  31-33  follows. 

The  main  effect  of  sediment  S  waves  for  small  S  wave  velocities  (100  m/.-ec) 
at  low  frequencies  is  to  add  a  significant  energy  loss  mechanism  through 
the  relatively  large  S  wave  absorption.  This  greatly  increases  the  RL 
and  results  in  the  large  peaks  in  Figs.  31  and  32.  At  higher  fre¬ 
quencies  (Fig.  33)  the  loss  to  sediment  P  wave  and  substrate  P  and  S 
wave  absorption  is  large  enough  to  produce  the  large  RL  in  the  absence 
of  sediment  S  waves  end  thus  the  effect  of  more  absorption  is  not 
important.  At  all  frequencies  the  effect  of  increasing  S  wave  velocity 
is  to  broaden  and  decrease  the  magnitude  of  the  interface  wave  peak 
until  it  merges  with  the  general  increase  in  RL  due  to  sediment  S  wave 
exc  l tat  ion. 
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A  more  realistic  layer  structure  would  include  a  shear  speed 
gradient.  The  additional  effects  of  a  constant  gradient  of  3  sec  1  is 
shown  in  Figs.  34-36.  This  gradient  approximates  that  expected  in  deep 
sea  marine  sediments.  The  surficial  shear  speed  was  taken  to  be 
100  m/sec,  giving  a  value  at  the  substrate  interface  of  400  m/sec. 

For  comparison,  RL  for  the  fluid  sediment  and  the  solid  sediment  with 
v^  “  400  m/sec  and  no  shear  speed  gradient  are  also  shown. 

As  seen  in  Figs.  34-36,  the  addition  of  an  S  wave  velocity  gradient 
appears  to  make  no  qualitative  change  in  the  effect  of  sediment  rigidity 
on  energy  lost  to  interface  waves.  The  shape  and  magnitude  of  the  RL 
for  a  homogeneous  layer  (dotted  curve)  fairly  well  approximates  that  of 
the  layer  with  a  gradient.  This  is  expected  since  both  sediment  structures 
have  the  same  S  wave  speed  at  the  substrate  interface.  The  major  effect 
appears  to  be  a  shift  in  the  location  of  the  P  wave  oscillation  in  RL. 

Sediment  parameters  typical  of  sands,  important  in  shallow  water 

continental  shelf  environments,  have  not  been  specifically  treated  in 

the  present  paper.  Sands  differ  from  clays  in  two  acoustically  important 

ways.  First,  the  sound  velocity  in  sands  is  higher  than  that  in  water 

(1780  m/sec  at  the  sediment  surface  compared  to  1540  m/sec  in  water). 

This  means  that  the  P  wave  is  evanescent  in  sands.  Specular  reflection 

from  the  water  interface  is  more  efficient  than  in  clays.  This  combines 

with  the  exponentially  decaying  evanescent  wave  to  determine  the  P  wave 

amplitude  at  the  substrate  interface.  Secondly,  the  gradients  of 

parameters  in  the  first  few  meters  are  substantially  larger  than  in 

clays.  The  P  wave  velocity  gradient,  for  example,  can  be  about  20  sec  ^ 

-1  22 

in  sands  compared  to  2  sec  in  clays.  These  large  near-surface 
gradients  raise  the  possibility  that  gradient  driven  P-S  coupling  may 
be  an  important  mechanism  for  exciting  S  waves  in  sands  at  low  frequencies. 
This  process  is  explicitly  neglected  in  our  computational  model. ^  If, 
in  fact,  this  additional  excitation  mechanism  is  negligible,  our  analysis 
procedure  could  be  extended  to  treat  sand  sediment  types. 


FIGURE  34 

REFLECTION  LOSS  versus  GRAZING  ANGLE 
FOR  A  100  m  THICK  CLAY  LAYER  AT  25  Hz 

THE  SEDIMENT  SHEAR  SPEED  AT  THE 
SUBSTRATE  IS  400  m/sec 


ARL:UT 
AS-79-1842-P 
PJV-  GA 
8-  17-79 


REFLECTION  LOSS  -  dB 


FIGURE  35 

REFLECTION  LOSS  versus  GRAZING  ANGLE 
FOR  A  100  m  THICK  CLAY  LAYER  AT  50  Hz 

THE  SEDIMENT  SHEAR  SPEED  AT  THE 
SUBSTRATE  IS  400  m/sec 


ARL:UT 
AS-79-1843-P 
PJV-  GA 
8- 17-79 
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- FLUID  SEDIMENT 

. SHEAR  SPEED  GRADIENT  =  0  sec-1 


FIGURE  36 

REFLECTION  LOSS  versus  GRAZING  ANGLE 
FOR  A  100  m  THICK  CLAY  LAYER  AT  100  Hz 

THE  SEDIMENT  SHEAR  SPEED  AT  THE 
SUBSTRATE  IS  400  m/tec 


ARL:UT 
AS-79-1844-P 
PJV-  GA 
8-  17-  79 
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V .  SUMMARY 


A  study  of  the  effect  of  sediment  rigidity  on  bottom  reflection 
loss  (RL)  from  typical  deep  sea  sediment  types  shows  that  sediment 
shear  (S)  wave  excitation  is  important  for  thin  sediment  layers  but  is 
negligible  for  thick  layers.  The  major  mechanism  for  S  wave  excita¬ 
tion  is  compressional  (P)  wave  conversion  at  the  sediment-substrate 
interface.  Little  energy  is  coupled  into  S  waves  at  the  water-sediment 
interface. 

The  mechanism  for  S  wave  excitation  provides  a  means  for  quantifying 
the  categories  of  thick  and  thin  layers.  It  is  the  P  wave  amplitude  at 
the  substrate  interface  that  is  important.  If  the  P  wave  amplitude  at 
the  substrate  interface  is  significant,  the  layer  is  twiels  if  it  is 
small,  the  layer  is  In  general  there  is  a  grazing  angle  0q 

separating  a  low  angle  region,  in  which  the  P  wave  turning  point  is  wel] 
above  the  substrate,  from  a  high  angle  region,  in  which  the  P  wave 
strikes  the  substrate  interface,  or  has  a  turning  point  near  the  interface. 
For  angles  less  than  0q  the  layer  is  thick  and  S  wave  effects  are 
negligible.  For  angles  greater  than  0q  the  layer  is  thin 
and  S  wave  effects  produce  significant  increases  in  RL.  For  physically 
thick  layers,  is  large  and  S  wave  effects  are  negligible  for  all 
grazing  angles,  except  for  a  shift  in  the  peak  structure  between  the 
critical  angles  for  substrate  S  and  P  wave  propagation.  For  physically 
thin  layers,  0q  can  be  0°  and  sediment  rigidity  will  be  particularly 
important  at  low  grazing  angles. 

For  thin  layers  the  dependence  of  RL  on  frequency  shows  the 
existence  of  high  and  low  frequency  regimes.  The  total  attenuation  of 
the  S  wave  traveling  one  way  through  the  sediment  provides  a  means  for 
separating  these  frequency  regimes.  For  total  S  wave  attenuation  greater 
than  20  dB,  an  S  wave  generated  at  the  substrate  interface  is  effectively 


absorbed  within  the  sediment  without  striking  the  other  interface.  Since 
we  assume  attenuation  is  proportional  to  frequency,  this  defines  a  "high" 
frequency  regime  in  which  the  propagation  of  sediment  S  waves  can  be 
neglected;  however  their  excitation  is  still  an  important  loss  mechanism. 

At  "low"  frequencies,  S  wave  attenuation  is  less  than  20  dB  across  the 
layer  and  the  additional  interference  effects  due  to  the  propagation  of 
S  waves  through  the  sediment  are  important.  Very  large  peaks  in  RL 
(>20  dB)  are  possible. 

A  theoretical  treatment  of  the  reflection  and  transmission  coefficients 
at  the  water  and  substrate  interfaces  by  means  of  an  expansion  in  the  small 
sediment  S  wave  velocity  provides  a  basis  for  developing  a  ray  picture  of 
the  effects  of  sediment  rigidity  on  RL.  The  basic  mechanism  for  S  wave 
excitation  emerges  as  P  wave  conversion  at  the  substrate  interface.  The 
dependence  on  S  wave  velocity  snd  attenuation  and  P  wave  attenuation  are 
predicted  by  the  ray  picture.  Trie.  interference  structure  is  found  to  be 
due  to  the  interference  of  a  wave  similar  to  that  generated  in  a  fluid 
sediment,  the  A^  wave,  and  waves  whose  amplitudes  depend  on  the  S  wave 
velocity  to  the  first  power,  the  A^  wave,  and  to  the  second  power,  the  A ^ 
wave.  The  parameter  dependencies  predicted  by  the  ray  picture  are  veri¬ 
fied  in  detail  by  results  obtained  from  the  computational  model. 

Parameter  studies,  with  ranges  of  variation  restricted  by  geophysical 
data,  are  used  to  identify  significant  parameters  in  determining  RL  in 
solid  sediments.  The  results  obtained  agree,  in  general,  with  the  results 
expected  from  the  ray  picture.  Important  sediment  properties  are  found 
to  be  density,  P  and  S  wave  velocities,  and  the  gradients  of  P  and  S  wave 
velocities.  These  must  be  specified  to  better  than  half  their  geophysi¬ 
cally  allowable  range  in  order  to  obtain  RL  to  within  1  dB.  Just  as 
important,  some  negligible  parameters  were  also  identified.  These  include 
the  gradients  of  sediment  density  P  and  S  wave  velocity  gradients,  all 
attenuations  except  the  sediment  S  wave  attenuation,  and  the  substrate 
density.  These  parameters  may  take  on  any  value  within  their  geophysical 
range  without  changing  the  reflection  coefficient  by  10%.  The  remaining 
parameters  can  be  important  under  some  conditions  and  negligible  in  others. 


Also  presented  were  the  results  of  parameter  studies  which  indicate 
the  accuracy  required  in  specifying  geoacoustic  parameters.  Parameters 
typical  of  a  deep  sea  sediment  layer  were  varied  by  ±10%.  The  resulting 
changes  in  RL  were  used  to  identify  parameter  sensitivity.  RL  from 
solid  sediments  was  found  to  be  particularly  sensitive  to  water  sound 
speed,  sediment  thickness,  density,  and  compress ional  wave  speed,  and 
to  substrate  shear  velocity.  Very  low  sensitivity  parameters  are  all 
attenuations  and  their  gradients,  and  the  gradients  of  sediment  density 
and  P  wave  velocity.  Other  parameters  produced  fractional  changes  in 
the  reflection  coefficient  of  the  same  size  as  the  fractional  change  in 
the  parameter. 

Sediment  rigidity  was  also  found  to  influence  the  energy  lost  to 
interface  waves.  The  Stoneley  wave,  generated  in  the  fluid  sediment 
case,  merges  into  generally  increased  RL  in  solid  sediments,  and  the 
peak  in  RL  at  a  particular  grazing  angle,  characteristic  of  interface 
waves,  cannot  be  identified  in  solid  sediments  with  realistic  geophysical 
parameters. 

In  summary,  a  computational  model  of  the  bottom  reflection  loss 
of  a  solid  sediment  layer  was  used  to  study  the  effect  of  sediment 
rigidity  on  bottom  reflection  loss.  For  thin  sediment  layers  compressional 
wave  conversion  at  the  substrate  interface  is  the  dominant  mechanism  for 
the  excitation  of  sediment  shear  waves.  Significant  increases  in  bottom 
reflection  loss  are  possible  at  fairly  high  frequencies  (hundreds  of 
hertz).  For  low  frequencies  the  effects  of  sediment  rigidity  can 
dominate  the  reflection  loss.  For  thick  sediment  layers  the  compressional 
wave  conversion  cannot  take  place  at  the  substrate  interface  because  the 
wave  is  refracted  relatively  high  in  the  sediment.  In  this  case  the 
effects  of  sediment  rigidity  are  negligible.  A  ray  picture  of  the 
acoustic  processes  was  developed  and  used  to  understand  and  interpret  the 
results  of  the  computational  model.  Parameter  studies  are  reported  which 
identify  important  bottom  parameters  in  determining  the  reflection  loss 
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in  solid  sediments.  The  influence  of  parameter  accuracy  in  predicting 
bottom  reflection  loss  is  treated  using  parameter  studies.  For  realistic 
bottom  parameters  the  loss  of  energy  to  interface  waves  does  not  appear 
to  be  an  important,  distinct  process,  as  in  the  case  of  a  fluid 
sediment , 
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